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APPLICATION  OF  CFAST  TO  SHIPBOARD  FIRE  MODELING 
III.  GUIDELINES  FOR  USERS 


1.0  INTRODUCTION 


There  are  two  primary  areas  in  which  computer  fire  modeling  could  be  of  benefit  to  the  US 
Navy  The  requirements  imposed  by  these  applications  are  different  and,  in  some  ways 
complementary.  The  first  application,  ship  design,  requires  a  capability  for  performing  detailed 
simulations  in  a  complex  environment,  but  has  few  restrictions  on  the  time  required  to  obtain 
results  For  the  second  area,  real-time  prediction  of  fire  behavior,  the  model  must  be  able  to  run 
taster  than  real  time,  but  the  predictions  may  be  less  detailed. 

These  requirements  approximately  correspond  to  the  capabilities  of  the  two  different  classes  of 
fire  models:  field  models  and  zone  models.  Field  models  typically  divide  the  region  of  interest 
(domain)  into  hundreds,  thousands  or  even  millions  of  small  volumes  (cells),  the  dimensions  of 
which  are  typically  on  the  order  of  centimeters  or  smaller.  Temperature,  pressure  species 
concentrations  and  other  variables  are  calculated  for  each  cell  as  a  function  of  time  and  each 
variable  is  represented  as  a  time  variant  scalar  or  vector  field  (hence  the  name  field  model).  The 
results  of  these  calculations  can  be  very  detailed,  both  in  spatial  resolution  and  in  terms  of  the 
amount  of  information  available  for  each  cell.  However,  field  models  are  very  slow  —  it  is  not 
unusual  to  require  hours  of  supercomputer  time  to  simulate  fractions  of  a  second  of  real  time 
even  for  physically  small  systems. 

At  the  other  extreme,  zone  models  commonly  use  only  one  or  two  very  large  homogeneous 
vo  umes  (zones)  to  represent  a  volume  that  might  require  hundreds  of  thousands  of  cells  in  a 
field  model.  In  addition,  zone  models  typically  ignore  physical  and  chemical  details  so  that  each 
zone  is  normally  represented  by  many  fewer  variables  than  would  be  the  case  in  a  field  model 
Hnally  zone  models  usually  represent  only  the  time  variations  of  the  variables  and  use  ordinary 
rather  than  partial,  differential  equations.  This  combination  of  factors  (fewer  computational 
elements,  reduced  numbers  of  variables  per  element  and  simplified  equations  for  each  variable) 
permits  zone  models  to  be  many  orders  of  magnitude  faster  than  field  models. 

Of  course,  this  gam  is  not  free  —  the  trade-offs  for  increased  speed  are  reduced  spatial  resolution 
and  lower  accuracy  of  the  predictions.  For  some  applications,  including  real-time  prediction  over 
a  relatively  short  time  span  (perhaps  30  minutes),  this  is  an  acceptable  compromise  and  it  is 
reasonable  to  use  a  zone  model.  In  other  applications,  the  lower  accuracy  and  resolution  are  not 
acceptable  and  field  models  must  be  used.  For  ship  design,  it  is  likely  that  there  will  be  a  role  for 

zone  models,  especially  in  the  early  phases,  but  that  field  models  will  be  needed  for  the  detailed 
design  phase. 


The  Consolidated  hire  Growth  and  Smoke  Transport  (CFAST)  model  (1]  is  an  example  of  the 
zone  model  class.  Developed  over  a  period  of  many  years  by  the  National  Institute  of  Standards 
and  Technology  (NIST),  CFAST  was  designed  to  provide  fast  and  reasonably  accurate 
predictions  ol  lire  and  smoke  spread  in  buildings.  In  recent  years,  the  US  Navv  has  funded 
improvements  in  the  CFAST  model  to  make  it  more  applicable  to  shipboard  fires.  With  Navv 
sponsorship,  CFAST  has  gained  capabilities  for  modeling  phenomena  that  are  absent  from  or  of 
little  significance  to,  building  fires.  This  has  included  mass  transport  through  vertical  vents 
(representing  hatches  and  scuttles)  [2],  energy  transport  via  conduction  through  decks  |3|  and 
improvement  to  the  radiation  transport  submodel  |4|.  Work  is  currently  underway  at  the  Naval 
Kcseaich  Laboialory  (NRL)  to  validate  the  most  recent  addition,  horizontal  compartmenl-to- 
compartment  heat  conduction  through  bulkheads. 

Duiing  the  development  ol  CbASl .  certain  simplifying  assumptions  (in  addition  to  those  which 
are  inherent  in  zone  models)  were  made  regarding  the  types  of  phenomena  and  the  complexity  of 
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important  because  that  was  the  least  complex  compartment  (it  was  a  rectangular  parallelepiped) 
and  allowed  us  to  simulate  the  fire  without  the  confounding  factor  of  geometric  complexity.  The 
fact  that  the  ventilation  ducts  were  sealed  significantly  reduced  the  "difficulty  of  the  problem, 
since  we  only  had  to  contend  with  buoyancy  flows  through  hatches  and  doors. 

The  work  proceeded  in  two  phases.  In  the  first,  we  developed  the  fire  specification  and  tested  it 
against  Laundry  Room  data  [6].  During  the  second  phase,  the  original  model  was  extended  by 
sequentially  adding  surrounding  compartments  while  keeping  the  fire  specification  unchanged 
[7],  In  order  to  minimize  bias  due  to  foreknowledge  of  the  actual  test  outcomes,  the  test  results 
were  not  used  in  either  phase  until  after  the  modeling  was  completed5,  at  which  time  the  model 
predictions  were  compared  with  test  results  to  determine  the  accuracy  of  the  model.  CFAST 
version  3.1 ,42,  was  used  for  this  project. 

The  purpose  of  this  report  is  to  summarize  the  findings  from  both  phases  and  to  present  them  in  a 
form  which  will  be  useful  to  engineers  or  designers  interested  in  applving  CFAST  to  shipboard 
fire  problems. 

2.0  CFAST  BASICS 

In  this  report,  we  discuss,  in  the  form  of  a  case  study,  the  application  of  CFAST  to  problems  of 
Navy  interest.  In  particular,  we  focus  on  the  problems  that  were  encountered  and  the  solutions 
which  were  developed  during  modeling  of  a  portion  of  the  ex-USS  SHAD  WELL.  We  start  with 
a  discussion  of  the  basics  of  using  CFAST  but  make  no  attempt  to  replicate  the  material  in  the 
definitive  CFAST  references,  which  are  the  NIST  user  guide  [8]  and  technical  reference  [9].  The 
former  discusses  CFAST  from  the  perspective  of  someone  who  simply  wants  to  use  CFAST  as  a 
tool  while  the  latter  provides  an  explanation  of  the  underlying  methodology  and  algorithms  and 
is  primarily  of  scientific  interest. 

2.1  The  CFAST  Input  File 

The  execution  of  a  CFAST  model  is  controlled  by  a  text  input  file  which  contains  all  of  the 
information  required  to  describe  the  desired  simulation.  Each  piece  of  information  appears  as  a 
parameter  associated  with  a  keyword.  Collectively,  the  keywords  and  corresponding  parameters 
constitute  the  vocabulary  of  the  CFAST  command  language.  This  implies  the  existence  of 
anothei  class  ol  limitations,  in  addition  to  those  previously  mentioned,  because  no  problem  can 
be  simulated  unless  it  can  be  described  using  the  CFAST  vocabulary.  W'e  will  encounter  several 
instances  which  illustrate  this  problem. 

Some  keywords  have  only  a  single  associated  parameter,  others  have  several  parameters  and  still 
others  may  have  a  variable  number.  Except  for  comment  lines,  which  begin  with  the  pound  sign 
(#)  and  are  ignored  by  CFAST,  a  keyword  appears  at  the  beginning  of  each  line  in  the  command 
file.  Keywords  may  be  included  in  any  order  and,  where  appropriate,  the  parameters  associated 
with  a  keyword  appearing  later  in  the  file  will  replace  values  that  appeared  earlier. 

For  the  following  discussion,  we  divide  the  keywords  into  four  categories:  ( 1 )  simulation  control; 
(2)  ambient  environment;  (3)  model  geometry  and  (4)  fire  description  The  keywords  which  we 


1  There  were  two  exceptions  to  this  rule.  First,  the  measured  fuel  mass  loss  rate  was  given  as  an  input  to  CFAST  to 
ensure  that  we  simulated  the  specific  fire  that  was  present  in  test  4_1()  and.  second,  we  used  an  experimental  value 
lor  the  CFAST  carbon  monoxide  production  parameter.  The  latter  was  done  onlx  after  it  was  determined  that  the 
carbon  monoxide  parameter  had  no  significant  effect  on  temperature  predictions  Futther  details  max  be  found  in 
rc  lore  nee  1 6  ] . 

Subsequent  to  the  date  ol  this  work.  C  I  AS  I  Version  4  was  released  lor  public  ihc.  The  major  new  feature  of 
Version  4,  horizontal  heat  conduction,  has  not  yet  been  validated  against  experimental  data 
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taforZti™  Wound'S  re“[8] APPCndiX  A  bU‘  **  is  a"  “^stive  list.  Further 

2.1.1  Simulation  control  keywords 

!"  "'hidl  «“  simulation  is 

oAheUr  S  the' T^rY  **  f“ 

mechanism  by  whicSTS^SSCSo? “  ‘ 


model .HI 

Listing  1.  Examples  of  the  CFAST  Simulation  Control  Keywords. 

f • • -«  provides  an  option 
on-screen  printouts  and  the  ime™, f  mle™'  ^ 
two  parameters  are  the  interval  between  cmnh  nnHoto  r  ^  , \  ’  c  n  seconds.  The  last 
copies  of  each  graph  that  should  be  produced^rMAder"  s“r?l'ds)  a,,d  lhe  number  of  hard 
_and  DUMPR  species  the  nrT^gSjl  '  “l  Wl’ich  evenls  °“ur 


fhe tnpu,rfi *  ^XruTtoZTins  im  „nru  ChalaC,er  r°"0wi"8  ll“  VERSN  keyword,  of 

cfaspt  2.x  wiii  ^r'"- 

to  specify  a  descriptive  title  (up  to  50  characters)  for  ihc  simulalii,  ERSK  ’  pCnni,s  ",c  uscr 

dim  ™lf  wiaraSSSn^^-  f'T"0"  ,ime'  is  «“  "umber  of  seconds 
depends  on  the  speed  of  the  comoute^md  thl  «L  ? ™n.nm"  time  of  the  model.  The  latter 
parameters  allows  the  user  to  set  the  number  r.(V  •  P  ^X,‘  ,  ie  scenari°-  The  second  and  third 

history  file  updates,  mspectheh  The  ^vn  i  ""“'f  ^  sccon,ds  be‘"  een  screen  updates  and 
screen  graphic  display  updates  and  the  desireH^n.  anietf,s  conliol  the  bine  interval  between  on- 

Rather  than  use  the  limited  built-in  graphics  capabilities  of  CFASt".'!^0/  CaCh  graphic  disPla>'- 
history  file  and  plot  the  results  off-line.  P  CFAST,  most  users  post-process  the 

explanation.  CFAST  allows  'the  userm  define' events0"^ ?  kcywords  and  squires  further 
control  line  the  erowth  of  the  lire  Events  e-m  n--  ”  ’  SUC  1  as.  opening  and  closing  vents  or 

lime  base  that  rs  used  by  ah  events^ defines  •USC,‘SLCCI  ied  !'nios  but  there  is  a  single 
since  the  start  of  the  simulation)  the  time  for  each  event  b>  IS,"lg  (,n  scconds  claPscd 

H*~  2.  for  vents  A  end 

;r' ‘A 

Ul,K  n,,Ks  l-  ,hrou?h  L  ™ng  the  FTIMF  keyword  and  then 
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specify  each  of  the  vent  openings  (as  a  fraction  of  the  maximum  width3)  for  each  of  these  times, 
even  if  there  is  no  change  for  that  vent  at  that  time.  This  vent  sequence  is  illustrated  in  Table  1. 


FTIME 

tl 

t2 

t3 

CVENT  A 

0 

1 

0 

0 

CVENT  B 

1 

1 

1 

0 

Table  1.  Example  of  the  CFAST  Event  Mechanism 

Events  can  be  triggered  at  the  elapsed  times  (in  seconds)  specified  by  FT1ME.  In 
this  case,  the  CVENT  commands  cause  vents  A  and  B  to  full v  "open  (1)  or 
completely  close  (0)  at  limes  zero,  t,,  t2  and  tv  Note  that  an  event  time  at  zero 
seconds  is  implied. _ 

There  is  an  implicit  t(l  in  FT1ME  but  events  must  be  explicitly  provided  for  all  times,  including  t„ 
(f.e.,  the  first  event  in  each  set  is  assumed  to  occur  at  t„,  even  though  zero  is  not  listed  in  the  time 
line).  This  is  because  CFAST  reads  the  events  from  left  to  right;  the  first  event  is  assumed  to 
occur  at  zero  time  and  subsequent  events  are  sequentially  associated  with  the  times  listed  in  the 
FT1ME  line.  If  an  event  is  missing,  then  the  following  events  will  be  associated  with  the  wrong 
time.  "  r 

Because  CFAS 1  uses  Imcui  intcipolation  to  calculate  values  lor  intermediate  times,  transitions 
occur  over  the  period  between  sequential  event  limes  rather  than  instantaneously.  If  the  transition 
is  too  steep  (/>.,  there  is  a  large  change  in  the  parameter  value  between  two  closely  spaced  lime 
points).  CFAST  may  encounter  numerical  problems  and  drastically  slow  down  or  even  stall.  If 
problems  of  this  nature  are  encountered,  it  is  recommended  that  the  transition  be  spread  over  a 
longer  period  of  time. 


(  I  AS]  only  permits  vent  widths  to  he  varied  This  is  appropriate  for  doors  but  does  not  correct l\  represent 
vertically  opening  vents  such  as  windows. 
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different  times.  ValueTtoXeT^  for 

SSn  ' ?o  JS  V?  ‘°  ““ 

switching  from  the  parameter  set  associated  with  the  firs,  simulate  to  thteffc  second"8'3""' 

example!' some  *of  the" hn^L  togMcPteK^beS' wite  h8®  *”  ‘he  ‘2pe  °f  fue!  mvolved  For 
instantaneously,  to  a  diesel  sprav  fire  This  was  m  h  i  a  uX^?e  Pan  ^ire  anc*  switched,  virtually 
parameters  characterist  P  o  Se  up  Jo  he  toe  oftS  nil  s,mu,a,inS  the  pan  fire,  using 
the  model  with  diesel  fuel  fire  parameKm  S  of  1 the  sPra>'  fire  then  restartinl 

2.1.2  Ambient  environment  keywords 

TAMB  indicates  the^n^  EAMR0?^  inhListing  2 

conditions.  There  can  onlv  be  one  instance  nfMch  !  !  egl0n’  EAMB  refers  to  the  exterior 

single  TAMB  line  implies  ^  only  a 

prior  to  the  start  of  the  simulation.  described  by  a  single  initial  condition 


#  Temp . 
TAMB  285.900 
EAMB  286.300 


Press . 
101300. 
101300. 


Elev. 
0.000000 
0 . 000000 


Listing  2.  Examples  of  the  CFAST  Ambient  Condition  Keywords. 

Site  r,;"3;  ‘“ssrrisr a,,d  w.-*™ 

kl22I!g3lla^^  mean  sea  level)  a,  temncmmmsanSSs  w^sS' 

above  mean'sea  level')3 ure'^i^SIi  ^leSevnUotroTtli'eTe r>aSCalS‘>  !'nd  relereilce  elevation  (meters 
can  correct  for  altitude-related  temperature  and  nil  ei(:Ilce  Pol™ls  needed  so  that  the  model 
elevation  can  be  important for hXr^  clianges-  Correction  for  differences  in 

shipboard  fire  modeling.  The  reference  elevation'can^eY1”1  expected  to  be  significant  for 
point  must  be  used  for  both  keywords  The  elev  uion  of  .l^i  Y  conven,ent  Point>  but  the  same 
and,  for  ships,  this  would  norma be  zero <ta  ievei  °  '°WCSI  COmPar",K"'  «  Wically  used 

X- W*  wind  speed,  the  elevation  a, 
These  parameters  are  only  significant  if  there  -ire  i,;,!|C'1  -CUJ‘UI|0n  °  vwnd  sPecds  at  other  heights, 
other  vents.  WIND  was  not  used  in  this  work-1  i  ^  1  Winds  1  iat  lnteracl  with  doors,  windows  or 
problems.  ^  th,s  wo,k  and  ls  not  «I«cied  to  be  important  for  most  naval 

2.1.3  Implicit  Geometric  Limitations 

complicated  geo.net  r!es.  IndhsTecfioiT' we’Siijscuv  rCS,ricfti?ns  0,1  our  “hiliiy  to  define 
following;  section,  tu  will  discuss  some  o|  these  restrictions,  including  the 


a-  al1  compartment  boundaries  are  rectangular; 
b.  each  compartment  lias  only  a  single  bulkhead,  which  wr 


ips  around  all  four  sides; 
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c.  for  any  compartment,  the  overhead,  deck  and  bulkheads  are  each  limited  to  a 
single  set  of  thermophysical  properties; 

d.  vent4  locations  are  undefined5; 

e.  all  horizontal  vents  are  rectangular;  and 

f.  vertical  vents  must  be  either  circular  or  square. 

The  first  limitation  is  due  to  the  fact  that  the  dimensions  of  each  compartment  are  described  by 
only  three  parameters:  DEPTH,  WIDTH  and  HEIGH.  This  implies  that  each  compartment  is  a 
rectangular  parallelepiped,  i.e.,  there  are  exactly  six  bounding  surfaces,  each  of  which  is  a 
rectangle.  The  deck  and  overhead  are  assumed  to  have  identical  dimensions  and,  since  there  is  no 
provision  for  defining  individual  walls,  CFAST  treats  the  compartment  vertical’ boundary  as  one 
continuous  entity6.  As  a  corollary,  there  is  no  way  to  specify  in  which  bulkhead  a  horizontal  vent 
is  located. 


Also,  each  boundary  has  only  one  associated  entry  in  the  thermophvsical  properties  database. 
Thus,  the  deck,  overhead  and  the  entire  wrap-around  bulkhead  each  has  only  one  set  of 
properties  —  it  is  not  possible  to  exactly  represent  compartment  boundaries  that  have  patches 
composed  of  different  materials7. 

These  inherent  limitations  force  us  to  make  approximations  in  order  to  describe  many  common 
shipboard  situations.  For  example: 

a.  many  compartments  do  not  have  simple,  rectangular  cross  sections; 

b.  standard  water  tight  doors  are  not  rectangles;  and 

c.  watertight  hatches  (except  for  scuttles)  are  usually  neither  square  nor  circular. 

Later  in  this  report  we  will  encounter  several  instances  where  approximations  were  required  to 
circumvent  these  limitations.  Examples  of  our  approximation  methods  will  be  presented,  along 
with  discussions  of  the  implications  of  these  methods. 

2.1.4  Geometry  keywords 

Examples  of  the  use  of  the  most  important  geometry  keywords  (see  Table  A-3)  are  given  in 
Listing  3.  Compartment  floor  elevations  (meters),  relative  to  the  reference  point  defined  by  the 
ambient  environment  inputs,  are  given  by  HI/F  and  the  dimensions  of  each  compartment 
(meleis)  are  specified  using  DEPTH,  WID1  H  and  HEIGH.  1  he  dimensions  of  each  compartment 


4  In  CPA  ST  terminology,  a  vent  is  any  opening  between  compartments  (or  between  a  compartment  and  the  outside), 
except  for  ventilation  duets,  which  have  their  own  special  set  of  keywords.  Since  we  did  not  use  ventilation  ducts  in 
these  simulations,  references  to  vents  appearing  in  this  report  refer  to  doors,  windows,  hatches,  scuttles  and  similar 
openings.  An  oddity  of  CFAST  nomenclature  is  that  vents  are  described  by  the  direction  of  the  flow  through  the 
vent,  not  by  the  orientation  of  the  vent  itself.  For  example,  a  door  is  a  horizontal  vent,  because  it  allow  s  horizontal 
How.  although  the  orientation  of  the  door  is  vertical.  Hatches  and  scuttles  permit  \enical  (low  and  are  classified  as 
vortical  vents. 

'S  I*oi  an  exception  to  this  rule,  see  the  discussion  of  the  HALL  keyword  in  section  2.1  4 

l  °r  «-*aeh  time  step.  CFAST  divides  the  wall  into  upper  and  lower  portions  and  independently  calculates  values  for 
the  stale  ot'  each  portion. 


It  is  possible  to  specify  boundaries  composed  of  multiple  layers  where  each  laser  may  have  different  properties. 
For  example,  a  bulkhead  could  be  specified  as  having  a  steel  core  with  cork  on  one  side  and  fiberglass  on  the  other. 
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-rdinaee  system  with  the  „rigi„  ,ocated  al  lhe 


#  Cmpt.  1  Cmpt .  2  Cmpt .  3 

#  Laundry  Psgwy  Wardrm 
#Floor  elevation 


Cmpt.  4 
NER 


Cmpt . 
CR 


Cmpt.  6 
Sail  1 


HI/F  0.00  0.00  2  1 

DEPTH  1.75  10.26  4*( 

#Y  dimen . 

WIDTH  6.07  2.22  8.! 

#Z  dimen. 

HEIGH  2.57  2.57  2.^ 

#Materials 

CEILI  SHIP3/8  SHIP3/8  SHI 
WALLS  SHIPLR  SHIPLRP  SHI 
FLOOR  SHIP3/8  SHIP3/8  SHI 
# 

# Laundry- Passageway  door 
#  Cmpt#  Cmpt#  Vent# 


Cmpt .  7 
Sail_2 

10.18 

0.91 


SHIP3/8  SHIP7/8 
SHIPLRP  SHIPWR 
SHIP3/8  SHIP3/8 


SHIP3/8  SHIP3/8 
SHIPNER  SHIPCR 


SHIP7/8 


SHIP3/8  SHIP3/8 
SHIPCR  SHIP3/8 
SHIP7/8  SHIP3/8 


SHIP3/8 

SHIP3/8 

SHIP3/8 


HVENT  1  2  1  o.6 

#  Cmpt#  Cmpt#  Vent#  Widt 

7ENT  1  2  i  i.oo 

#Passaceway-Wardroom  hatch 

#  Cmpt#  Cmpt#  Area  Type  n 

WENT  2  3  0.78  ^ 

# Laundry -Wardroom  heat  conduction 

#  Cmpt#  Cmpt# 

CFCON  1  3 


Width 

0.66 

Width@t0 

1.00 


Soffit 

1.90 

Width@tl 

1.00 

circular ; 


square) 


Listing  3.  Examples  of  the  CFAST  Geometry  Keywords. 

These  are  the  major  keywords  used  to  define  the  geometry  for  a  tvniral  pfa qt  a  ,  • 

is  the  compartment  floor  heiohl  ('relative  m  ,kA.V  y  typical  c  FAST  simulation.  HI/F 

WIDTH  and  HEIGH  specify The  dime  sfo ns  C„  fference.elev^*on  from  Listing  2);  DEPTH, 
identify  the  ceiling,  S’  Ld  F  %  COnipartmeills-  CELI  WALLS  and  FLOOR 

contained  in  a  separate  database  file  THERMAL  DF  H V FN T  c° P ° ' r ° S  ,°f  -,hcSC  ,11aterials  are 

WALLS  and  FLOOR, "i^clivd  b>  "'c  ke',"ords  CEIL1, 

construction  materials  as  they  appear  in  the  m  iteri  .1  °P  1CSC  keywords  are  lllc  names  of  the 
default,  the  file  THERMAL  DF  is  used  \  ^‘tabase  which  is  a  text  file.  By 

to  specify  a  different  rUename  ' The mop  J  ’1"  k  **  b>'  usi"g  the  THRMF  keyword 

list  of  thermophysical  ponies "oX  c h^  n  ,he  "a™  of  lhe  material  and  a 

thickness  and  emissivity.  The  term  "thermophysica  properties^ Xlme"i  Spcc,,,c  hcal-. density, 
because  the  values  apply  to  a  specific  lhiekness  V  ,  l  7  ,11  ls  somewhat  misleading 

S  “* includc 

"'ickncss)  can  he  defined  for  ,„c 

ddlerent  materials  (or  different  thicknesses  of  the  sat, tc  nnterhlLtV,"!',  T  '°SI<’"S  c"l"""clcl1 

s.une  malutal),  then  some  approximation  must 
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be  applied.  In  a  subsequent  section  of  this  report,  we  discuss  some  methods  for  making  these 
approximations. 

It  is  important  to  emphasize  that  CFAST  treats  walls  as  wrap-around  entities  so  that  each 
compartment  is  surrounded  by  a  single  wall,  not  by  four  separate  walls.  Thus,  not  only  is  it 
impossible  to  exactly  model  situations  where  a  single  wall  has  regions  of  different  materials, 
some  approximations  must  also  be  made  in  those  cases  where  different  walls  have  different 
compositions. 

Each  of  the  keywords  mentioned  above  (except  for  THRMF)  is  followed  by  a  list  of  values,  one 
value  per  compartment.  Compartment  numbers  are  determined  by  the  order  of  the  parameters, 
i.e.,  the  first  parameter  following  the  keyword  applies  to  compartment  one,  the  next  to 
compartment  two  and  so  on.  There  is  no  keyword  to  declare  the  number  of  compartments; 
instead,  CFAST  determines  this  value  by  counting  the  parameters  associated  with  these 
keywords.  If  N  parameters  are  found,  CFAST  assumes  there  are  N  compartments.  Furthermore, 
it  also  assumes  that  "compartment"  (N+l)  refers  to  the  world  outside  of  the  model  domain. 

It  is  important  that  the  same  number  of  parameters  appear  in  each  of  these  lines.  If  different 
numbers  appear  on  different  lines  of  the  input  file,  CFAST  will  run  without  reporting  any  errors 
but  the  results,  as  might  be  expected,  are  likely  to  be  erroneous. 

HVENT  and  WENT  are  used  to  define  horizontal  and  vertical  vents,  respectively.  HVENT 
requires  six  parameters;  in  addition,  there  are  three,  seldom  used,  optional  parameters.  The  first 
two  required  parameters  identify  the  two  compartments  that  are  connected  by  the  vent.  As 
mentioned  above,  compartment  numbers  are  determined  by  the  order  in  which  their  dimensions 
were  declared.  Recall  that,  for  an  N-compartment  simulation,  compartment  (N+l)  represents  the 
ambient  environment.  This  convention  allows  the  definition  of  vents  that  connect  between  a  real 
compartment  and  the  external  environment. 

CFAST  permits  up  to  four  HVENTs  between  the  same  pair  of  compartments,  so  the  third 
HVENT  parameter  is  a  number  (from  one  to  four)  that  identifies  which  particular  vent  is  beinu 
defined.  The  next  three  parameters  specify  the  vent  width,  the  height  of  the  soffit  (top)  and  the 
height  ol  sill  (bottom),  respectively. 

The  first  optional  parameter  is  the  cosine  of  the  angle  between  the  wind  vector  and  the  vent 
opening8.  This  is  applicable  only  to  vents  that  connect  to  the  exterior  and  is  meaninsful  only  if  a 
wind  speed  is  declared  via  the  WIND  keyword.  The  last  two  parameters  define  the  horizontal 
position  of  the  vent  with  respect  to  the  origins  of  the  two  connected  compartments.  These 
parameters,  which  were  recently  added  to  CFAST  to  permit  modeling  of  flows  in  long,  narrow 
corridors  using  the  HALL  keyword,  will  not  be  discussed  further  because  they  were  not  used  in 
our  work.  Additional  information  about  the  HALL  keyword  may  be  found  in  reference  1 10], 
which  reports  on  the  validation  of  the  new  corridor  flow  algorithm. 


Like  HVENT,  WENT  also  requires  two  parameters  to  define  the  connected  compartments.  In 
addition,  the  area  of  the  vent  must  be  given  and  the  shape  is  specified  as  either  circular  or  square. 
Vents  of  other  shapes  must  be  replaced  by  equivalent1'  round  or  square  openings.  This  is  an 
example  of  a  case  in  which  the  CFAST  vocabulary”  restricts  the*  user’s  ability  to  define  the 
problem. 


8  the  orientation  ol  lire*  vein  is  defined  by  the  outward  lacing  normal  vector. 

’’  In  this  context.  "equivalent"  refers  to  the  How  resistance  of  the  vent  and  not  Hr  the  actual  vent  area.  However,  for 
rents  that  do  not  deviate  too  much  from  circular  or  square,  adjusting  the  vent  area  is  probably  sufficient. 
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CVENTES  lSf0ntf  7“  c'osure  events'  The  parameters  for 

HVENT)  and  a  list  of  fractional  vent  openings  The fiMfarfa'10i"  "umber  (as  defined  by 
zero  time  and  subsequent  entries  correspond m t  JiL  •  fract,on  refers  to  the  vent  opening  at 
CVENT  was  included  in  the  SHADWELL/6K8  <^ces®lY®  eveat  times,  as  established  by  FTIME 
always  set  to  1 .0  because,  in  the parameter  was 
open  during  the  fire.  There  is  no  corresponding  Svu  ^Hf  moddm^  a11  of  the  vents  were  fully 
possible  to  define  closure  events  for  vertical  vents  Y  VertlCal  VentS  S0’  at  Present> *  is  not 

'»  fHe  vertical  beat  conduction 

doefnlt SS  “X  of  the 

be^an  important  factor  in  most  muhi-deck,  shipboardTemS  dufto 

was  performed  butlSVa"' k“  not  avai,^,e  at  Ihe  time  this  work 
Validation  of  this  new  feature  has  not  yet  been  nerfnL  j”0! ‘°  allow  for  this  possibility, 
horizontal  conduction  algorithm  is  noTknown  perf°rmed’  50  the  accuracy  and  utility  0f  the 

2.1.5  Why  specify  the  fire? 

p  ,hre  iOUS  W"v 

fuel  and  allow  the  model  to  calculate  the  behavior  of  the  file  a!Vde«!ops?  ‘yPC  a"d  omoun'  of 

sx":  rciplcr is  **  <«  ■<« 

reactions  and  any  aliempt  10  model  them  would  requir/deiSrf''  7  01  ln(eraclln8  chemical 

thermodynamics  and  kinetics  of  those  reactions  Mn  ?f!  L  d  ?  1 formation  regarding  the 

a  field  model,  which  would  defeat  the  purposes  of  CFAST^"161113  y*  "  Would  rccluire  the  use  of 

the  fire.  In  faciJUs  somewhat  no  attemPl(°  actually  calculate 

ejects  of  the  fire  but  does  not  actually  mSiel  the  firelSf  R  ,7  T**  CFAST  modc,s  ^ 
chemical  reactions,  CFAST  uses  user-supplied  1  :  f*  R'  ?thei r  than  a«empt  to  simulate 

concentrations  of  various  combustion  products  for  •.  o  'S  ,al-,os  to  es"ma(e  the  relative 

of  fuel  consumed  is  based  on  the  user-specified  fire  siz^descripliom  °‘  fucL  The  amounl 

CFaS-^S1  imTmt  -  d-  -  of 

specification  and  considerable  effort  must  be  exnenffm  !,!  °  Para[11etcrs  "lvolvod  in  (he  fire 

consistent.  Furthermore,  if  the  goal  of  the  model  is  m  l  T'0  lhat  the  description  is  internally 

~  1,1  Ul:  W°rk)'  ^  Is  requSi  lo  lUl 

is  the  most  important  singlVpmXnoVofihTllI/^  rc,casc  rale  of  (lie  fire 

dial  (lie  heal  release  rale  will  be  the  most en  ' 1  Thercr°,v-  ^  ^P«:t 

^ -S2S  «’a= 

S  ^o'auii'ihri^^  r' 
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CFAST  provides  a  capability  for  specifying  "object"  fires  which  are  intended  to  represent  real 
objects,  such  as  furniture,  that  may  ignite  during  the  course  of  the  scenario.  Depending  on  the 
user  inputs,  object  fires  ignite  at  a  specified  time,  surface  temperature  or  surface  heat  flux.  Object 
fires  were  not  used  in  our  work  and  will  not  be  discussed  further. 

2.1.6  Fire  specification  keywords 


The  most  common  fire  specification  keywords  are  given  in  Table  A-4  and  Listing  4  provides 
examples  of  their  use.  The  location  of  the  main  fire  is  given  by  FPOS  and  LFBO.  FPOS  specifies 
the  coordinates  of  the  fire  within  the  compartment,  using  the  compartment  coordinate  system  that 
was  previously  described.  LFBO  specifies  the  number  of  the  compartment  in  which  the  fire  is 
located. 


#  X  Y  Z 

FPOS  C . 91  1.83  0.19 

#Fire  Cmpt 

LFBO  1 

#Fire  Type  (1  =  unconstrained;  2  =  constrained) 

LFBT  2 

#Mass  pyrolysis  rate 

FMASS  0.0253  0.0229 

#  Kol  Wt  Rel  Hum  LOL  He  Init  T  Ign .  T  Rad.  tract. 

CHEMI  184.  100.  10.  4.19E+007  285.9  330.  0.30 

#H:C  mass  ratio  (fuel  composition) 

HCR  0.143  0.143 

#0 : C  mass  ratio  (fuel  composition) 

02  0.0  0.0 

#HCN:fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HCl:fuel  mass  ratio  (pyrolysis) 

HCL  0.0  0.0 

#Toxics : fuel  mass  ratio  (pyrolysis) 

CT  0.0  0.0 

# Soot : CC2  mass  ratio  (combustion) 

OD  0.06  0.06 

#CO : C02  mass  ratio  (combustion) 

CO  0.056  0.056 

CJET  OFF 

Listing  4.  Examples  of  the  CFAST  Fire  Specification  Keywords. 

The  fire  is  located,  at  the  coordinates  given  in  FPOS,  in  compartment  one  and  is  a  constrained 
(type  2)  fire.  The  fuel  mass  loss  rate  starts  at  0.0253  kg/sec  and  declines  to  0.0229  kg/sec  at  1250 
seconds.  CHEMI  specifies  the  fuel  molecular  weight  and  heat  of  combustion,  as  well  as  several 
other  parameters.  HCR  and  02  specify  the  amount  of  hydrogen  and  oxidizing  agent  in  the  fuel, 
respectively;  OD,  CO,  HCN,  HCL  and  Cl  define  the  production  rates  lor  various  pyrolysis  and 
combustion  products.  CJET  OFF  invokes  the  standard  ceiling  jet  convection  algorithms. _ 


Fires  may  be  defined  as  either  unconstrained  (type  one)  or  constrained  (type  two)  by  using  the 
LPBr  keyword.  These  types  control  the  way  in  which  the  fire  will  respond  to  reductions  in 
atmospheric  oxygen  concentration.  In  an  unconstrained  fire,  CFAST  always  uses  the  specified 
heat  release  rate;  for  a  constrained  fire,  the  heat  release  rale  will  he  reduced  below  the  user- 
specified  level  if  there  is  insufficient  oxygen  to  support  that  combustion  rale.  We  should  also 
note  that,  for  the  unconstrained  fire.  CFAST  does  not  calculate  any  species  concentrations.  The 


permit  uncons,rai"ed  ™ 
S"'  ™ASS  Sd  CHEMLThe  ptrj^ass^ 

time-dependent  -  that  ,s,  they  correspond  to  specific  times  defined  by  FT1ME  °  316 


dQ/di  =  aHc  dm/di 


Eqn.  1 


the  hea,  ^omb^TOe  ?a« ert  'the  Canine  a"d  ^  iS 

Obviously,  if  all  three  parameters  are  given,  the  fire  is  over-SDecified  pfa  qt  ™i 

of  the  three  and,  if  all  three  are  present  uses  nniv  toot  ?  sPdCined-  Ob  AST  only  requires  two 

three  keywords  in  the  order  FQDOT  FMASS  and  CHFMTrFA^?301]1?  6’ lf  ?  USer  includes  the 
and  4H„  to  caicnlate  new  FQDOT  val’uefin^accordance^vilh  ^afionT"  "*  ™ASS 

m„,addrn  ,0,therhena'  °f  «"»•»»**«».  other 
humiditv.  the  lower  oxygen  limit  fLOIi  of 7hl f i' ecu,ar  we,ght  of  the  tuel,  the  initial  relative 

vapor  ignition  tempeSe'S  ft?  ridiLwe  tZt  “mSceTT^  T,™ 

parameters  requires  further  explanation.  ’  significance  of  some  of  these 

monoxid^arl'rela^vefy “  0X»“-  diraide  *"d  carbon 

water  vapor  However, 

concentration  to  use.  The  relative  humiditv  ic  ,iwh  .  .11  n?eds  f°  be  told  what  initial 

modified  by  the  water  contribution  from  romhnct-  nff.  CU  ale  1 1IS  initia|  v:due,  which  is  then 
infrared  radiation,  its  concentration  affects  rid  h  fin  1  ecausc  wale,‘  vaPor  is  a  strong  absorber  of 
addition,  humidity  is  a  factor  in  ,he 

wMcIm'ik  fire1 ^iil'beUfinvmS1' TOsfimhf"  7’CCn!nUi?"  <in  ™'ume  percent)  below 
(constrained)  fire  Tlte  sunprclsio,  fnncTi  „  k  c"lorced  0nly  for  the  case  <>f  a  type  two 
LOL  +  0.0 1  down  to  zero  extinguishment)  at  lheVor°Tl’  Va.'yrln®  ,ro'."  onc  <no  suppression)  at 
is  a  compromise,  made  by  I Ik  CFAST  dove '7 °  n Tlle/ deraidl  '  due  lor  LOUs  10%.  This 
typically  reported  to  he  irf  f he  h  T"  1“  ,hal  L0L  “ 

approach  zero  in  some  lai”c-sc  ile  fires  M?i  p  ,  ,  c  .  es  hut  has  been  reported  to 

turbulence,  which  cou”d  cSe  oxtlen  ied  nL  Ls  ?Sl'Tb  y’  ,h'.S  discrePa,lc>'  «  related  to 
die  nominal  oxygen  concentration  i.x'too  low  to  support  bm'n'ing';0'” '  ’US"°"  cvc"  "'h,'n 

\"C  firC  P,UT  or-  by  radiation,  hot,,  the 
is  important  because  i,  amS?  *"  P^uced.  the  radiative  traction, 

and  ceiling.  The  default  value  k  "oo  . ,  m  s  i  '  lcmPC';"urcs  ol  the  gas  layers,  the  walls.  n,n.r 
clean  fires  (tnctltanol.  for  example),  this  value's  uk!  high  CbSS  '*  °UCl)  filcs  ^  vcr-v 
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There  are  seven  other  keywords,  in  addition  to  CHEMI,  that  affect  the  chemistry  of  the 
combustion  process:  HCR,  02,  HCN,  HCL,  CT,  OD  and  CO10.  The  first  two  of  these  are  related 
to  the  fuel  composition,  the  next  three  to  the  behavior  of  the  fuel  during  pyrolysis  and  the  last 
two  to  the  combustion  process.  Collectively,  these  keywords  provide  the  information  used  by 
CFAST  to  calculate  the  production  of  various  pyrolysis  and  combustion  products;  the 
concentrations  of  these  products  are  then  tracked  as  they  spread  through  the  simulated  system 
The  parameters  associated  with  all  of  these  keywords  are  treated  as  events  so  it  is  possible  to 
specify  fires  in  which  the  effective  chemistry  changes  with  time. 

HCR  and  02  are  the  hydrogen-to-carbon  and  oxygen-to-carbon  mass  ratios  of  the  fuel, 
respectively.  HCR  is  used  to  determine  the  amount  of  water  vapor  produced  when  the  fuel  is 
burned  and  02  reflects  the  amount  of  oxidizer,  if  any,  that  is  incorporated  into  the  fuel.  Any 
oxygen  provided  within  the  fuel  itself  reduces  the  need  for  atmospheric  oxygen  and  therefore 
affects  the  heat  release  rate  and  oxygen  consumption  in  constrained  fires.  The  02  keyword 
applies  only  to  rocket  fuels  and  other  special  cases.  Fuels  which  contain  oxygen  within  the 
molecular  formula  may  be  considered  to  be  partially  pre-oxidized  (methanol  can  be  thought  of  as 
partially  oxidized  methane,  for  example).  In  these  cases,  the  effects  of  the  oxygen  are 
automatically  accounted  for  by  the  heat  of  combustion  and  the  02  keyword  is  not  needed. 

HCN  specifies  the  ratio  of  the  mass  of  hydrogen  cyanide  in  the  pyrolyzate  to  the  mass  of 
pyrolyzed  fuel.  Similarly,  HCL  defines  the  mass  ratio  of  hydrochloric  acid  to  fuel.  CT  is  similar 
to  HCN  and  HCL  in  that  it  controls  the  rate  of  production  of  a  pyrolysis  product  and  this  product 
is  tracked  as  it  spreads  through  the  system.  However,  CT  is  not  a  real  product,  rather,  it  is  a 
quantity  which  represents  (based  on  empirical  correlation  functions)  the  total  toxic  hazard  of  the 
actual  products.  It  takes  into  account  the  typical  toxic  effects  of  such  things  as  carbon  monoxide 
acid  gases  and  oxygen  depletion.  The  CT  input  affects  the  predicted  concentration  of  "total 
toxics"  but  does  not  affect  the  prediction  of  temperature  or  of  any  of  the  real  fire  products. 

The  last  two  keywords,  OD  and  CO,  specify  mass  ratios  for  products,  but  these  ratios  are 
calculated  with  respect  to  the  mass  of  carbon  dioxide  in  the  product,  rather  than  to  the  mass  of 
the  fuel.  OD  is  the  soot-to-carbon  dioxide  ratio  and  CO  is  the  carbon  monoxide-to-carbon 
dioxide  ratio.  OD  strongly  affects  the  predicted  concentration  of  soot,  of  course,  but  also  has  a 
major  indirect  effect  on  temperatures  because  soot  is  usually  the  dominant  term  in  calculating 
radiative  energy  transport. 

The  OD  and  CO  ratios  vary  with  combustion  conditions  and  the  user  is  responsible  for  providing 
physicalh  reasonable  values.  This  can  be  difficult  lor  two  reasons:  (1)  the  parameters  are  ratios 
of  products  and  any  change  in  the  combustion  process  that  affects  one  product  is  likely  to  affect 
the  other,  making  it  difficult  to  intuitively  estimate  the  net  effect  on  the  ratio;  and  (2)  due  to 
sensitivity  to  the  combustion  conditions,  it  is  often  not  possible  to  find  literature  values 
appropriate  to  the  simulation  conditions.  Furthermore,  even  if  values  consistent  with  the 
expected  conditions  are  specified,  the  conditions  calculated  by  the  model  may  be  verv  different 
from  those  expectations. 


For  example,  suppose  that  you  anticipate  that  the  fire  will  always  be  well  ventilated  and  specify 
an  OD  value  that  is  appropriate  for  this  case.  If  the  model  predicts  some  unanticipated  ventilation 
ellect  (perhaps  a  vent  closure  far  from  the  fire  has  a  much  greater  effect  than  was  anticipated), 
then  the  model  may  throttle  back  the  lire  (assuming  it  is  a  constrained  fire).  It  is  likely  that  the 
predicted  soot  production  will  then  be  too  low  (because  CFAST  will  not  increase  soot  production 
to  account  lor  the  reduced  oxygen)  and.  since  temperatures  are  strongly  affected  by  soot 
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In  common  usage.  ()1)  ictcis  lo  optical  density  while  C  ()  and  02  relcr  to  the  species  carbon  monoxide  and 
oxygen,  respectively.  (TAST  assigns  different  meanings  to  these  terms  and  this  non-standard  usage  sometimes 

ate  always  used  in  the  C  bASI  sense  and  species  are  spelled  out. 


causes  contusion.  In  this  report.  Ol),  CO  and  02 
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reduTet^hea^rdea^'ra^wil^be^o^loweMp8  may,be  Sr°SSLy  in  OTOr'  The  dire“  ««*  <* 

change  in  soo,  proDerti«7n  STT*!”  “A™™*  of  «>e  fire,  but  the 

increase,  decrease  or  remShe  sarnf  $°  ,hat'  fUr‘her  fr0m  ,he  fire’  '<™P«dures  could 

c„CRha",d  02>  th?se  parame,ers  are  l«ri«ic  to 
does  not  calculate  a  self-consistent  fire& iStU  °n  conditions.  However,  because  CFAST 

calculate  the  composition  of  the  pyrolysis  product^bas^on^h^e^user^nputs011  ^  ^  ^ 

transfer  to  the  cemng^uringTheTnitial  stagesof  afire  t STf"2  m°d-el  °f  conductive  heat 
"splash"  against  the  ceilin<*  spreading  radiallv  frnm  -’f  ?  ^ire  P  u]?e  nses  rapidly  and  may 
jet  can  hetTt  the  ceiling  more  efficientfv Than  h/X  P,°lnt  ab°Ve  the  fire'  This  s°-caIled  ceiling 
of  later  times.  This  effec?  can  be  trnnortam  Vf***'  gas that  is  more  ^ 

becomes  insignificant  as  the  ceiling  temperature  approaches  thaTo? th^oas^  Hflre,,,bu'  normall>' 
becomes  thick  enough  to  act  as  a  buffer  against  the  rising  fire  pTume  “  ‘he  Upper  layer 

disabiled^br  XEIL/NG"R(iea,coifvect^n*  enabledbf f  ",?FF "  fceilinS  j«  convection  algorithm 
CFAST.  the  values  "WALL"  aT“ALL"  were  -wf?"?g)'  In  some  earlier  versio"s  of 
convection  to  walls;  "ALL"  enabled  both  ceiling  fnd  ,5ctlvated  the  jet  algorithm  for 

wall  jet  algorithm  produced  erroneous  resnlt^in/f  2  JCtS'  H?wex'er’  11  was  f°und  that  the 
available  in  CFAST  3  14  The  CJETninmpiP  d,  consequently,  that  feature  is  no  longer 

the  model  is  the  same  as  if  "OFF  had  to  snenTiS'  f  T  f*  S  fTALL"  but  the  behaviorV 
as  the  "CEILING"  parameter  Spec,f,ed'  S,m,larly’  ALL  Produces  l|ie  same  effect 

were  nmuTcd  ^nlhls  wL*  SfTd  be  although  they 

position  specified  by  FPOS.  FPOS  is  a  constim  hn^Pmr’u'6  base  of,tlie  fIames  above  the 
keyword  provides  a  mechanism  for  varying  »|1P  »i  *1  treated  as  an  event,  so  this 

could  be  used  si muime  a  deS. r° ph" itl  ^1 ‘!oh  the  b“sc  of , h<f  n .T T ‘ F°r  “ 

down.  Since  our  simulations  used  pool  fires  in  which  ,iU  i  !  ie.f  e  decreases  as  the  pile  burns 
to  the  heigh,  of  the  eompartmem  tfc  keywo’rd  “needS  °'  "*  »“  ■*«*»*  -hi . . 

sophisticated  pluni'c ''niodeK'n'1  wddeh^rainment0^  "  "as  Provided  so  that 

incorporated  into  CFAST  However  those  nlumo  mnrt  |S  a  functlon  °J  fire  area>  could  be 
present.  FAREA  is  not  used  P  models  were  never  =>dded  to  CFAST  and,  at 


3.0 


development  of  a  fire  specification 


Submarine ^Ven, tab X  “TK  ""  specticatio.a  for  the  particular  case  of 
expect  to  be  folhtwed  in  apn  vin  .  CFAST  ,  '.S,pioc'ss  .«?  s"™lar  to  that  which  we  would 
factor:  in  order  to  co  , , are  o ur  msuhs  udth  ‘ hypothcfcal  new  ship  design,  except  for  one 
contras,,  the  ship  dcsi  "wr  L?no  a  “ Lr  “  ,di'  a’  T  l,ad  10  si"«'hne  a  specific  fire.  In 

representing  many  possible  casualties,  uouldlrc  sitmdmcd  dtuT^Z X'ign  pow  ™'C'y  "'CS' 


3.1 


Donning  the  Laundry  Room  Model 


development  of  the  fire  specification. 

can  he  run.  we  must  define  all  of  the  nomine  paramcX 
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the  fire  compartment.  For  this  case  study,  the  fire  compartment  was  the  Laundry  Room,  which  is 
illustrated  in  Figure  3.  The  simplicity  of  this  compartment  minimizes  the  difficulties  in  setting  up 
those  parameters  and,  in  fact,  that  was  one  of  the  reasons  for  selecting  test  4_10  as  the  subject  for 
this  case  study.  The  CFAST  input  file  describing  this  single-compartment  simulation  is  shown  in 
Listing  5. 


VERSN 

3  SHADWELL/688  Laundry  compartment  1 

# 

Sim . time 

Print  Hist.  Disp.  Copies 

TIMES 

1250 

13  0  0 

# 

to  tl 

FTIME 

1250  . 

DUMPR 

model . HI 

# 

Temp . 

Press.  Elev. 

TAMB 

285.900 

101300.  0.000000 

EAMB 

286.300 

101300.  0.000000 

a 

Cmpt .  1 

# 

Laundry 

# Floor 

elevation 

HI  /F 

0 . 00 

dimen. 

DEPTH 

1 .75 

#Y  dimen. 

WIDTH 

6 . 07 

#Z  dimen. 

HEIGH 

2 . 57 

# Mater ia Is 

CEILI 

SHIP3/8 

WALLS 

SHIPLR 

FLOOR 

SHIP3/8 

Lisl 

ing  5.  CFAS  r  Input  File  lor  the  Laundry  Room  Simulation 

This  simple,  one-compartment  model  of  the  Laundry  Room  was  used  for  development  of  the  fire 

spool  1 1C 

anon. 

15 


#  Cmpt# 
HVENT  1 

#  Cmpt# 

CVENT  1 
#Fire  Cmpt 
LFBO  1 

#  X 
FPOS  0.91 


Cmpt# 

2 

Cmpt# 

2 


Vent# 

1 

Vent# 

1 


Width 

0.66 

Width@t0 

1.00000 


Soffit  Sill 

1.90  0.00 

Width@tl 
1.00000 


He 

4 . 19E+007 


Y  Z 
1.83  0.19 

Type2(1  =  unconstrained;  2  =  constrained) 

#Mass  pyrolysis  rate 
FMASS  0.0253  0.0229 

Mol  Wt  Rel  Hum  LOL 

CHEMI  184.  100.  10. 

#H:C  mass  ratio  (fuel  composition) 

HCR  0.143  0.143 

#0 : C  mass  ratio  (fuel  composition) 

02  0.0  0.0 

#HCN : fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HC1 : fuel  mass  ratio  (pyrolysis) 

hcl  O.o  0.0 

#Toxics : fuel  mass  ratio  (pyrolysis 
CT  0.0  0.0 

#Soot:C02  mass  ratio  (combustion) 

00  0.06  0.06 

#C0:C02  mass  ratio  (combustion) 

C0  0.056  0.056 

CJET  OFF 


Init  T 
285.9 


Ign .  T 
330. 


Rad.  fract. 
0.30 


Hjling  5  (corn'd).  CFAST  Input  File  for  H,e  I  . . .  silmilalion 


3.1.1 


Simulation  control 


-  11“^  —  <!«  sinm  Union  lime  (TIMES, 

begun  ,o  die  somewhat  before ihl,  si ? teL™KZProX,",;"cly  '?;<’* *c0"ds'  bul  the  lire 
lime-relaled  factor  Ilia!  dirccllv  fceled  he  r  "  ,  “  ’  ’  wf  s,;l  10  i2'’°  seconds.  The  only 

lest.  As  is  discussed  below  his  cliani  L  ,s  f  '  r°  S'°W  chanSc  ">  mi!ss  loss  rale  during  ilie 

loss  rates  only  for  the  start  and  end  ofthe  test"  Accordin«!y  fS^!niCient  lo  fPccif>'  mass 

at  1250  seconds  (recall  that  there  is  nn  imnlir'ii  «  f-  ^  defines  a  single  event  time 

(DUMPR)  was  given  a  oeneric  name  MOI1FI  ri/6111.1"116  Ul  Ze,°  seconds)-  T,le  output  file 
simulation  began  at  ti.nefero  MODbL.HI,  and  no  restart  file  was  needed  since  these 

3.1.2  Ambient  environment 

uinhurair'tompc^Sfo?^  “!c  ««•  «  used  aclual  pre- 

sliipboard  tcmpcraliircs  do  not  vary  widely  from  comiviitoiei  i  m'-T"1  lempcralul>-s-  Typically. 
«l  specific  lempernture  info, rmaiiin.  i,  wlnllTbc  aTs  '  to  % T pT-'A f  “bSC"C<; 

h^^iuS;  ,^s',m':ucr  .la  *r fircs- 

s;,,,  1,1  any  casi.  I„c  Ulcc.s  ^  “  g  !^tol7ty 


for  the  test  ease,  ambient  temperatures  were  deter 
temperature  readings  during  the  period 


.'term mod  by  taking  the  average  of  all  air 
pnoi  to  ignition,  flic  values  for  the  Laundry  Room 
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(285.9  K)  and  the  Laundry  Passageway  (286.3  K)  were  used  for  the  initial  temperature 
parameters  for  TAMB  and  EAMB,  respectively.  The  rationale  for  using  the  passageway,  rather 
than  the  exterior,  temperature  for  EAMB  was  that,  from  the  perspective  of  the  fire  compartment, 
the  passageway  js  the  exterior.  That  is,  air  entering  through  the  door  will  be  at  passageway 
temperature  and  conduction  through  the  walls  will  be  governed  by  the  temperature  difference 
between  the  laundry  and  the  passageway. 

Of  course,  conduction  through  the  deck  and  overhead  areas  involves  other  compartments,  which 
may  be  at  different  temperatures.  Most  of  the  test  area  was  at  about  the  same  temperature  (within 
several  degrees),  so  this  is  not  likely  to  introduce  a  significant  error.  In  any  case,  the  limitation  of 
being  able  to  specify  only  a  single  external '  temperature  makes  this  type  of  approximation 
unavoidable. 

Pressures  for  both  TAMB  and  EAMB  were  set  to  101300  Pa  (one  atmosphere)  and  the  reference 
elevations  were  set  to  zero  (sea  level).  As  was  mentioned  previously,  the  absolute  elevation  is 
significant  only  when  modeling  high-rise  buildings.  For  a  model  in  which  the  maximum 
elevation  differences  are  only  on  the  order  of  a  few  tens  of  meters,  these  parameters  are  not 
important. 

3.1.3  Laundry  Room  geometry 

The  Laundry  Room  deck  elevation  was  set  to  zero  (HI/F),  which  implies  that  the  TAMB  and 
EAMB  reference  points  are  at  this  level.  Since  this  compartment  is  a  rectangular  parallelepiped 
we  can  immediately  specify  the  DEPTH,  WIDTH  and  HEIGH  parameters  as  1.75  m  (5.74  ft)’ 
6.07  m  (19.9  ft),  and  2.57  m  (8.43  ft),  respectively.  Note  that  we  have  adopted  the  convention 
that  width  (y-axis)  is  parallel  to  the  length  of  the  ship  and  that  depth  (x-axis)  is  in  the  athwartship 
direction.  Because  this  compartment  meets  the  CFAST  constraint  that  all  surfaces  are  assumed  to 
be  rectangles,  we  are  able  to  represent  the  geometry  exactly,  with  no  approximations. 


All  of  the  boundaries  of  the  Laundry  Room  were  constructed  of  steel  of  various  thickness.  Since 
the  SHAD  WELL/688  lest  area  was  heavily  sooted,  the  entries  for  steel  in  the  thermal  properties 
database.  THERMAL.DF,  were  edited  to  change  the  emissivily  from  the  default  value  (0.9), 
which  is  t\ pical  ol  clean  steel,  to  1.0.  which  is  a  reasonable  value  lor  steel  coated  with  carbon 
black. 

The  Laundry  Room  overhead  (CEILI)  and  deck  (FLOOR)  were  constructed  of  0.95  cm  (0.38  in  ) 
steel,  which  is  one  of  the  values  (SH1P3/8)  in  THERMAL.DF.  However,  the  bulkheads 
(WALLS)  were  made  of  several  steel  plates,  of  different  thicknesses,  welded  together.  Since 
CFAST  does  not  permit  the  use  of  multiple  materials  for  any  single  boundary,  we  defined  a 
fictitious  material,  SHIPLR,  to  represent  these  bulkheads.  Because  die  plates  were  made  of  the 
same  material,  we  used  the  standard  thermophysical  properties  (density,  conductivity  and  heat 
capacity,  for  example)  of  steel  and  varied  only  die  thickness.  To  ensure  that  the  effective  heat 
capacity  of  the  model  bulkhead  was  the  same  as  the  total  of  (he  actual  bulkheads,  w  e  required 
that  out  bulkhead  mass  equal  to  the  actual  mass.  I  his  implied  that  the  material  volumes  would 
also  be  equal  and  it  followed  that  the  area-weighted  mean  was  the  proper  average  to  use  for  the 
thickness.  Based  on  this,  SHIPLR  was  specified  to  have  a  thickness  of  0.0076  m  (0.30  in.).  The 
relevant  portions  of  the  modified  version  of  THERMAL.DF  are  shown  in  Table  2. 

The  joiner  door  leading  into  the  Laundry  Passagew-ay  was  defined  as  an  HVENT  eonneetins 
compartment  one  (Laundry  Room)  to  compartment  two  (the  exterior  of  the  modeling  domain^ 
This  w»as  vent  number  one  and  had  a  width  of  0.66  m  (2.2  ft),  a  soffit  of  1.90  m  (6.23  ft)  and  a 
sill  at  zero.  The  vent  was  always  fully  open,  so  the  parameters  of  the  corresponding  CVENT 
were  set  to  1 .0  for  both  event  times. 
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!S5S  S  S  ! 

SHIP5/8  48  550  £.0127  1 

t”L  4S  I  I  1  g  |  oj)Q76  |  1  [ 

Table  2.  Material  Property  Database  Entries  for  Steel  Plate. 

for  .he  Laundry  Room 

in  the  SHADWEIJL  test  area6 Th^SHtP^p  vanous  thicknesses  of  steel  plate  used 

simulation;  the  other  entries  are  standa^d^hicknerses  of  steel  %eClf,Cal.1X  for  this 

set  to  one  (rather  than  the  default  0  9)  because  ,SS1VIties  were 

covered  with  soot.  oecause  surfaces  in  the  test  area  were 

3.2  The  Fire  Specification 

m  Xl's^m  W  a99CftxP6.M  2!  ^  ^  *itf!  the  P“  “»>ered  a,  0.9,  m 

mass  loss  measurements,  elevated  the  nan  0  IQ  m  comPartl^eiit  origin.  The  load  cell,  used  for 
with  the  FPOS  keyword  and  LFBT  specified  i  m  £ec^-  These  coordinates  were  given 

oxygen  depletion  Luld  belnd^  ‘7>  S° ^  the  effects 

was  not  included  because  the  current  version  of  PRaTt/  °f  ,he  fire’  087  m2  (9-37  f{2) 

previous  discussion  of  the  FAREA  keyword).  AST  d  n0t  Use  t  lIS  in,ormation  (recall  the 

In  CFAST,  burning  is  treated  as  a  two-step  process. 

Fuel  — >  Fuel  vapor  +  Pyrolysis  products  j£^n  2., 

Fuel  vapor  +  Oxygen  -*  Combustion  products  Eqn  2b 

Where  Equations  2a  aud  2b  represent  pyrolysis  and  combus.ion,  respectively. 

as  rue,,  pyrolysis  or  combus.iou 
parameiers  are  considered  in  IhffoC^  Lee  see, ioT^  Cy  ^  These  '"ree  »f 

3.2.1  Fuel  parameters 

.bCCaUSC  Sr-d  only  on  ,he  nauue  o,  ,„e 
properties  or  mos,  real  riels  (as  Zd '„>  dwT?,P'°CCSu  OIISC"inS  ,his  “  ■'«  feel  that  ,he 
have  high  batch-to-batch  variability  and  few  Ivin-I,  ^i'  7  pUIt-  surr°ga,es-  such  as  hexane) 
Values  for  many  fuels  may  often  be  obmine  lm m  m"  n°.  ^aracten/.ed  in  any  significant  way 
military  and  standards  organization  specifications"  and  simihr^sf"0 ll,dm£  chcnucaI  handbooks, 
m  m‘  ,ha'  many  °f  ,hCSC  SOUrCCS'  staiukirds  dc^men^ 


Eqn.  2a 
Eqn.  2b 


I  he  American  Society  for  Testing  ami  Materials  (A STM)  is  one  example  of  such  a  standard 


s  organization. 
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limits.  As  a  result,  literature  values  must  be  considered  to  be  estimates,  rather  than  exact  values, 
for  the  properties  of  any  given  fuel  batch. 

Fortunately,  this  problem  is  usually  not  critically  important  because,  in  most  cases,  we  have  no 
way  of  knowing  the  details  of  some  hypothetical  future  fire  and  therefore  must  consider  a  range 
of  possibilities.  It  is  reasonable  to  start  with  nominal  values  taken  from  these  references,  on  the 
assumption  that  they  will  be  representative  of  typical  batches  of  the  given  fuel.  Simulations 
should  then  be  run  using  a  range  of  deviations  from  the  nominal  values  to  account  for  the 
possibility  that  some  fuels  will  be  atypical  —  out  of  specification  or  contaminated  with  other 
fuels,  for  example. 


In  our  case,  the  fuel  was  a  marine  diesel,  but  no  further  details  were  given.  We  used  values  taken 
from  Table  5.2  of  Kanury  [13]  for  the  molecular  weight  (184  gm/mol)  and  for  the  heat  of 
combustion  [4.19  x  107J/Kg  (1.80  x  104  BTU/lb)].  The  vapor  ignition  temperature  was  estimated 
from  the  minimum  flash  point  for  number  two  diesel  fuel,  which  is  given  as  52  °C  in  the  ASTM 
specification  for  diesel  fuels  [14].  Since  all  refiners  provide  a  safety  margin  above  this  minimum 
requirement,  we  assumed  that  the  actual  ignition  temperature  would  be  10%  higher,  or  57  °C 

(330  K).  All  three  of  these  are  specified  as  parameters  of  the  CHEMI  keyword. 

According  to  the  fuel  specification  [15],  the  minimum  hydrogen  content  of  F-76  diesel  fuel  is 
12.5%  by  weight.  Using  this  value,  and  the  molecular  weight  of  the  fuel,  we  estimated  HCR  to 
be  0.143.  We  ran  simulations  with  a  range  of  HCR  values,  as  discussed  in  detail  in  reference  [6], 
and  found  that  this  parameter  had  only  a  limited  impact  on  the  model  predictions  in  most  cases. 
These  fuel  characteristics,  and  the  sources  for  the  values,  are  summarized  in  Table  3.  Finally^ 
since  the  fuel  used  in  the  SHADWELL/688  tests  had  no  included  oxidizing  agent,  the  oxy«en- 
carbon  ratio  (02)  was  zero. 


Property 

Value 

Source 

Molecular  Weight  (gm/mol) 

184 

[11] 

Heat  of  Combustion  (J/Kg) 

4.19  x  10A7 

[11] 

Ignition  Temperature  (K) 

330 

[12] 

Min.  Hydrogen  Fraction  (%) 

12.5 

[13] 

I _ Table  3.  Typical  Properties  for  Marine  Diesel  Fuel. 

3.2.2  Pyrolysis  parameters 


The  pyrolysis  process,  illustrated  in  Equation  2a,  involves  evaporation  or  thermal  decomposition 
of  the  (assumed)  solid  or  liquid  fuel1-.  During  this  process,  some  fuel  components  evaporate  to 
foim  a  fuel  vapoi,  others  may  break  apart  into  lower  molecular  weight  species  which  then 
evaporate  as  more  fuel  vapor  (this  is  typical  of  some  plastics  and  other  polymers)  and  still  others 
decompose  directly  to  form  products  without  undergoing  combustion,  the  latter  mechanism 
typifies  the  production  of  acid  gases,  such  as  hydrochloric  acid  and  hydrogen  cyanide. 

Because  the  goal  ol  this  work  was  to  determine  how  well  CFAST  simulated  a  specific  test  fire,  it 
was  very  important  that  the  fire  description  accurately  represent  the  actual  fire.  Therefore,  the 
pyrolysis  rate  (l-'M  ASS)  inputs  were  calculated  from  actual  test  data.  To  reduce  the  experimental 


Some  luels  (methanol,  lor  example)  evaporate  to  form  a  flammable  vapor:  other  (such  as  plasties)  undemo  a 
chemical  decomposition  (pyrolysis)  process  first  so  that  the  vapor  is  chemically  dissimilar  to  the  original  fuel. 
CI-AST  does  not  distinguish  between  these  two  eases  —  I'M  ASS  represents  the  total  mass  lost  from  the  original 
luel.  without  regard  to  whether  thermal  degradation  was  involved. 
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time,  as  shown  in  Figure T We'found^hat  th^S  slld!Jjg,  ave,ra^®)  and  Pitted  as  a  function  of 
curve  shown  in  the  figure  (R2  =  0  955)  w  ch  w  i  ?  C‘°f fitted  with  *e  exponential 

The  latter  was  very  nearly  linear  (R2  >  0.9998)  an7since  CFAST^1116  T'5  l0$S  rate  Curve- 
between  event  times,  FMASS  was  defined  hv  tL  iff  CFAST  uses  linear  interpolation 

[0.0253  and  0.0229  kg/s  (0.0558  and  0.  oiosd| b/^.'^s^dviTy f. rateS  a'  “*°  and  1250 


/cm.  1  ML  SIXL,<-S-  consequently,  we  also  set  CT  equal  to 

II  1 1  u  oi  o  known,  or  siKnt'riiwl  ilnt  iim  in  o  ..  .  ,  .  ... 

(for  example,  if  it  had  been  a  poIvvinvIehloVidic'ii  ;s,^ni],canl  anuninls  ol‘  nitrogen  or  chlorine 

III  fiBd  ;t*~-  “  ^ 
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either  maximum  or  minimum  allowable  concentrations  (sometimes  both)  for  contaminant  species 
(e.g.,  nitrogen  compounds  in  fuels);  in  those  cases,  it  would  be  reasonable  to  estimate  values  for 
HCN  and  HCL  following  a  procedure  similar  to  that  used  for  HCR. 

3.2,3  Combustion  parameters 

The  combustion  process,  shown  in  Equation  2b,  is  usually  the  most  difficult  to  handle  because 
the  parameters  involved  are  a  function  of  the  combustion  conditions,  which  we  do  not  know  a 
priori.  The  primary  parameters  in  this  category  are  the  production  factors  for  soot  and  for  carbon 
monoxide  (OD  and  CO,  respectively). 

CFAST  calculates  the  reduction  in  combustion  due  to  oxygen  depletion13  but  does  not  account 
for  the  changes  in  the  combustion  product  ratios  that  also  occur  when  the  fire  is  oxygen  starved. 
Thus,  we  are  in  the  position  of  trying  to  guess  what  the  burning  conditions  will  be  so  that  we  can 
provide  the  correct  product  ratio  inputs,  when  those  very  inputs  influence  the  burning  conditions. 
As  previously  mentioned,  the  inability  of  CFAST  to  dynamically  adjust  the  combustion 
parameters  as  a  function  of  instantaneous  conditions  is  probably  the  most  significant 
shortcoming  of  the  model.  Unfortunately,  accurate  combustion  calculations  would  require  the 
use  of  a  field  model  incorporating  detailed  reaction  kinetics  and  that  w  ould  negate  the  advantages 
of  a  fast,  simple  zone  model. 


1  *  This  assumes  that  a  constrained  lire  lias  been  specified,  as  was  the  case  in  our  work. 
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Temperatures  were'founlT tolVTighl^  °D’  ranginS  from  zero  to  0.10. 

relatively  insensitive  to  the  exact  value  nfnnf  1  °n  c  Presence  or  absence  of  soot  but 
chose  0.06  as  an  abouI «•  <*■  Ultimately  we 

(methanol,  for  example),  it  would  be  important  to  use  f value  neajzerm  CXpeCted  t0  burn  cleanly 

a"d  literature  is  of  lit,, e 

expected  to  be  about  0.013  [16]  The  CFAST  default  v^f  o  n'?  rnonoxlde  to  carbon  dioxide  is 
value  for  a  "typical"  fire  which  would  Fndudl  non tamS S 3  reasonable  ">«" 
we  determined  that  a  CO  parameter  value  of  0.056  was  SaWe  foAesI  TJT  ^  5)> 

parameter  is  not  criticaho^he  fire^pec^cation  and'iTma^be^ac^^t'h1?2  in  man>’  cases.  this 
near  0.05.  However,  if  the  carbon  monoxide  concentradon^rf  P  5  C  t°,!,se  an  arbitrary  value 
does  become  critical.  Without  a  correct  value  for  thifoaramet^rn  ’/hen  the  C0  Parameter 
be  very  inaccurate.  At  present,  there  does  not  aonear  ,conlcentratlon  predictions  will 

best  approach  is  probably  to  run  the  simnlatirm  t0  be  a  g(lod  soIutlon  t0  this  problem.  The 
the  results  as  a  ra?,ge  of  pVs^Se  leentmi"  ns  *  Vme»  °f ' differe“  CO  inpuS  and  presem 

3.2.4  Miscellaneous  parameters 

fall  into  any  of  the  aboTe^al^  calculations  but  do  not 

humidity,  the  lower  oxygen  limit"  the  ^  Son  and  teaEr's^T*™'  ** 

temperature,  on  the" ass ulnpl i on 'tha^ the 1  fuXpans  were^  ^  '°  thei  laUMdry  room  ?mbicnt 

with  the  compartment  prior  to  ignition  of  (he  fuel  pproximately  in  thermal  equilibrium 

.ha, 

Default  \ alues  lor  (he  lower  oxv^pn  limit  a  m  \  j  . . 

respectively)  were  used  because  ilhose  wcre  belieielno  h‘p ‘al'Ve  fra?tion, inpuls  (30<*  and  10%, 
no  better  information  was  available  The  radiation  fr  ip?  b  reasonable  values  for  diesel  fuels  and 
the  fire  plume  and,  i„  real  fires !  that  is' cmurolled  °,n  Jhe  .a™u.'«  of  soot  in 

saluiated  fuels  typically  produce  cleaner  fires  ih-m  h  i  ^  i  /  16  nature  °*  lbe  luel  (for  example, 
amount  of  available  oxygen.  ‘  d  highly  unsaturated  fuels)  and  partly  by  the 

^^h  CF.A.ST ,  the  amount  of  soot  is  controllcvl  hv  oo  * 

radiative  fraction  should  be  consistent  with  OD  The  <\  »f  n^r  ’  *!S  Pre V1°usly  discussed,  and  the 
for  "dirty"  fuels  but,  for  "clean"  fuels  sud^i  “S  te1,or‘hfse  Paranieiers  are  appropriate 
should  be  reduced.  ‘  ethanol),  both  the  radiative  fraction  and  OD 

simulation.  Since  w^wea^  prinvln'l v' i 'n 'Vlic f  * °C ^  few  seconds  of  the 

(fens  of  minutes),  we  set  CJET  to  "OFF"  for  (his  work  m°(  °  p,edlctK)ns  ,or  long  time  periods 


Results  of  the  Fire  Specification  Development 


The  Laundry  Room  was  modeled  (usina  the  fire  sneeifie  uinn  . . .  ,  ,  , 

. . . . . . -  .n-<>Upio  r t;,- 


here  is  to  illus,rate  the  general  performa"“ ,hat  ™v  ■* 

However,  prior  to  making  this  comparison,  it  was  necessary  to  determine  which  thermocouDles 
measured  the  upper  layer  and  which  the  lower  layer.' In  the  following  paragraphs  we  briefly 
discuss  the  experimental  data  processing  methodology  and  explain  how  the  thermocouples  were 
apportioned  between  the  layers.  F 

The  Laundry  Room  was  instrumented  with  a  single  thermocouple  tree,  illustrated  in  Figure  6 
that  had  thermocouples  at  elevations  of  0.14  m  (0.46  ft),  0.50  m  (1.64  ft),  0.95  m  (3.12  ft)  1  27 

m  u  ’  17  *'^,m  (6-37  ft)  and  2.50  m  (8.20  ft)  above  the  deck.  The  temperature  reported  by 

each  thermocouple  was  recorded  at  one  second  intervals  for  several  minutes  prior  to  ignition  of 
the  test  fire  and  the  Laundry  Room  pre-ignition  temperature  was  calculated  as  the  ensemble 
average  (the  mean  of  all  measurements  for  all  thermocouples).  Zero  point  offset  errors  were 
corrected  by  adjusting  the  readings  so  that  the  mean  pre-ignition  temperatures  for  each 
thermocouple  was  equal  to  the  ensemble  mean. 


The  Cb AS  I  -predicted  interface  height  (Figure  7)  was  then  used  as  the  demarcation  between 
owe i  and  upper  layer  thermocouples.  Ignoring  the  transients  that  occurred  during  the  first 
several  seconds  the  interlace  was  found  to  lie  at  approximately  0.6  m.  Accordingly  the  bottom 
two  thermocouples  were  assigned  to  the  lower  layer  and  the  lop  four  to  the  upper  layer. 


For  each  layer,  the  appropriate  thermocouples  were  averaged  and  the 
deviation  was  estimated.  As  seen  in  Figure  8A  (upper  layerfand  f  igure 
predicted  tempci  alines  were  very  close  lo  the  measured  values  (well 
deviation )  except  at  the  start  of  the  fire.  During  that  early  period, 
overestimated  the  air  temperatures. 


population  standard 
SB  (lower  layer),  the 
within  one  standard 
CTAST  consistently 
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Figure  8A.  Comparison  of  Actual  and  Predicted  Upper  Layer 
Temperatures  lor  the  Laundry  Room 

The  error  bars  represent  one  standard  deviation  in  the  test  data.  The 
assignment  of  thermocouples  to  the  upper  or  lower  layer  was  based  on  the 
interface  height  estimated  by  CFAST. _ 
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Figure  8B.  Comparison  of  Actual  and  Predicted  Lower  Laver 
Tempeialures  for  ihe  Laundry  Room 

4.0  DEVELOPMENT  OF  A  GEOMETRY  SPECIFICATION 

except  for  1 1 1  c  c  orrcctlo  n‘ '  fo  ^  he  'b  u’l  kYil'a  d'l !  iVck  n  c  s  s  HoPlC  ‘hat  "°  aPP«;oxi  mat  ions  were  needed, 
required  significant  simplified ion  be  We  "!“'<*  the  olhcr  compartments 

compatible  with  CFAST  limitations  geometry  could  be  specified  in  i  manner 

^  "  be  significant:  I)  non- 
Laundry  Passageway  provided  the  bes  e^  mn  dmfb'fuo“? ^rmophysical  properties.  The 

rectangular  geometry  and  also  provided  in  ex-i?nb?f  W  d'  ICllltlcs  ,n  baling  with  no n- 

The  Wardroom  added  vertical  JXton  ,o  Z  ST"  "  Tf  a.verti«" 

we  encountered  a  case  in  which  the  standard  thermonhvli  ’  N lgal,on  hM»'pment  Room, 
this  section,  we  illustrate  the  approximation  hJ  ?  ,  '  Properties  were  not  applicable.  In 

Ihe  details  of  the  modeling  of  these  compartments  |,,v  |d  ‘"ld  Pr°NCnt  1,10  rationales  for  each, 
i  be  repeated  here.  compa.lmcnts  have  been  reported  previously  |7|  and  will 


4.1  Problems  of  Noil-Rectangular  Geometry 


As  may  be  readily  seen  from  Figure  9  the  deck  nl-m  «r  m; . 

s"ai,ci1  . . . ■  - . . 
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outboard  bulkhead  is  cut  away,  producing  a  five-sided  vertical  cross  section.  Any  one  of  these 
factors  would  make  it  impossible  to  create  an  exact  description  for  CFAST.  In  addition,  the 
bulkheads  were  made  of  several  different  thickness  of  steel  plate. 


There  w  ere  three  horizontal  vents  and  one  vertical  vent  in  the  Laundry  Passageway.  One  of  the 
horizontal  vents  was  already  included  in  the  model  as  part  of  the  Laundry  Room;  the  other  two 
were  standard  Navy  water  light  hatches  having  a  0.66  m  (2.17  ft)  w  idth  a  0.23  m  (0.75  ft)  sill 
and  a  2.04  m  (6.69  ft)  soflit.  Unlike  the  Laundry  Room  door,  these  hatches  had  rounded  corners, 
but  the  radius  of  curvature  was  negligible  (on  the  order  of  a  few  centimeters).  All  of  these  vents 
u'ere  fully  opened  during  the  tests,  so  each  of  the  fractional  openings  (defined  by  the  CVENT 
parameters)  was  set  to  1.0.  The  new  vents  were  defined  as  connecting  compartments  two 
(Laundry  Room)  and  three  (exterior).  For  consistency  with  the  Laundry  Room,  the  deck 
elevation  was  set  to  zero. 

In  this  section,  we  present  some  approaches  to  approximating  compartment  dimensions.  Issues 
related  to  vertical  transport  and  thermophysical  properties  will  be  addressed  in  subsequent 
sections.  How'ever,  before  considering  the  details  of  our  methods,  we  will  discuss  some 
considerations  that  may  be  more  broadly  applicable. 

4.1.1  General  considerations 

We  previously  noted  that  CFAST  was  designed  with  the  assumptions  that  each  compartment  has 
exactly  six  bounding  surface  and  that  each  surface  is  rectangular.  These  assumptions  are  implicit 
in  the  equations  used  to  calculate  surface  areas. 

A,llVk=  A„,hil  =  DEPTH  *  WIDTH  Eqn.  3a 
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Eqn.  3b 


^ lower  blkhd  2  *  (DEPTH  +  WIDTH)  *  I 
Aupperbikhd=  2  *  (DEPTH  +  WIDTH)  *  (HEIGH  - 1)  Eqn  3c 

oftlie  compwtment.1  °f  ‘he  be‘Ween  the  '°Wer  a"d  zones'  ™^sured  from  the  deck 

Likewise  the  equations  for  the  zone  volumes  are 


V lower”  DEPTH  *  WIDTH  *  I 
Vuppcr  =  DEPTH  *  WIDTH  *  (HEIGH  - 1) 


Eqn.  4a 
Eqn.  4b 


selfZsten,  Howeve,  a  prob  fZ  ?"elep!ped’  3  *"d  4  are  correct  and 

In  that  case,  three  dimenltion  pa  ameters  are  no°  sXienTfo  »n2 rcfla"8ular  Parallelepiped, 
inputs  can  be  adjusted  so  th/some 

inputs  that  can  make  all  of  them  simultaneously  mrme.  ua  1 th  I  ,  combination  of 
to  satisfy  and  which  to  sacrifice?  y  eCL  S°'  hoW  do  "e  declde  wl’ich  equations 

sprcification^wliidi  defui^^he'producti^n'rate'ffo  ’l,ere  ^  larSel>’  delermi"ed  by  ">e  fire 
volume.  The  latter  dictates^row  tnuch^liliuion^occtirs'aiid.^hereforeiaffecls^he  conce^f^l^0111 

T  and  transport  play  imporunfrofe"5 

example)  and  a  rettfal,^^  temf  rature  or  concentration,  for 

area,  respecttvely.  Thus,  a"d 

m  s  senses  tifz  ho"\ ,he  rvr 

Accordingly) vie  im^id Jil>rnaidered-  "  ‘°  be  ciuica'Hfacwr! 

chosen  so' as  to  give  app, oxtmate  compartment  dimensions  be 

because  they  directly  influence  only  the  trattspoti^uadons  CTC  ’ed  ‘°  bC  °'  'eSS  imporlance 

smmmsmm 

IlIlllP“===f=g 


Note  that  vertical  veins  ilo  not  have  this  desire 
the  upper  layer  of  one  compartment  and  the  km 


v  of  complexity.  Because  they  al 
e*  htyet  ot  another,  thev  can  neve 


u 

'I 


a>s  reside  at  the  interlace  between 
span  layers. 
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Unfortunately,  given  the  constraint  of  a  fixed  volume,  surface  areas  and  compartment  heights  are 
not  independent  parameters.  By  fixing  volume  and  height ,  we  have  implicitly  fixed  the  overhead 
and  deck  areas  as 


Adeck  =  Aovhd  =  Vacl  /  Hw  =  DEPTH  *  WIDTH  Eqn.  5 

where  Vact  and  Hac[  are  the  actual  volume  and  height  of  the  compartment.  The  only  remaining 
option  is  to  specify  either  DEPTH  or  WIDTH,  the  other  being  determined  by  Equation  5. 

This  approach  works  well  if  our  primary  interest  is  in  horizontal  fire  spread  and  there  are 
horizontal  vents.  However,  if  the  primary  goal  is  to  simulate  vertical  conduction  through  decks, 
then  it  is  important  that  the  overhead  and  deck  areas  be  accurate.  In  that  case,  a  better  procedure 
might  be  to  conserve  compartment  volumes  and  overhead  areas.  The  compartment  height  would 
then  be 


HEIGH  =  Vacl/Aaclovhd  Eqn.  6 

where  Aact  ovhd  is  the  actual  area  of  the  compartment  overhead.  This  again  provides  some 
freedom  in  selecting  DEPTH  or  WIDTH  but  we  must  realize  that  horizontal  vent  flows  will  be 
less  accurate  with  this  approach. 

In  either  case,  the  choice  of  DEPTH  and  WIDTH  is  somewhat  arbitrary.  In  our  work,  we  used 
the  actual  value  of  one  of  these  parameters,  if  an  unambiguous  value  was  available.  Since  the 
version  of  CFAST  available  for  this  project  does  not  support  horizontal  heat  conduction  between 
compartments,  we  did  not  make  a  special  effort  to  correctly  specify  bulkhead  areas.  We  did, 
however,  compare  the  actual  and  approximated  surface  areas  in  order  to  estimate  the  amount  of 
error  introduced  into  the  model. 

Because  vents  directly  influence  both  mass  and  energy  transport,  correctly  specifying  the  vent 
parameters  was  also  important.  Vent  areas  and  shapes  were  chosen  to  be  as  accurate  as  possible. 
In  our  simulations,  vertical  vents  always  required  some  approximation  since  the  hatches  were 
rectangular  rather  than  square. 

4.1,2  Estimation  of  actual  volumes 

In  determining  the  Laundry  Passageway  dimensions,  our  first  step  was  to  estimate,  as  accurately 
as  possible,  the  actual  volume  of  the  compartment.  Due  to  the  complexity  of  the  Laundry 
Passageway,  this  process  was  somewhat  complicated. 

As  shown  in  Figure  10,  the  Laundry  Passageway  was  divided  into  two  parts.  Section  1  and 
Section  2,  such  that  the  actual  compartment  volume  was  conserved 


V  =  V,  +  V,  Eqn.  7 

where  V  is  the  total  Laundry  Passageway  volume  and  V,  and  V,  are  the  volumes  of  Section  1 
and  Section  2,  respectively. 
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Section 

1 


Figure  10.  Conceptual  Division  of  the  Laundry  Room  Passageway 
I - into  Sections _ 

For  Section  1,  the  volume  is 


V,  =  H,  *  W,  *  D, 

and  the  dimensions  may  be  read  directly  from  Figure  9 

H  i  =  2.57  m 
W,  =  2.44  in 
D,=  1.75  m 

giving  a  volume  of 


lie]  n.  8 


Lqn.  9a 
Eqn.  9b 
Lqn.  9c 


1  '  /  in 


1 1,  as  a  function  of  the  location  alon o  'Thc'lV i dt l/ ax i , J’cu LiT^Th ^  \>,K| '  area’  sh‘nvn  in  Fi"urc 
cross-section  to  get  "  *  calculate  the  \oIuine,  we  integrated  this 


Vj  =  Ja(w)  d\\>  =  47.50  111'' 
where  the  width  of  Section  2,  taken  from  Figure  9.  is 

W,  =  8.51  m 


Lqn.  I  I 


Lqn.  1 2 
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Figure  11.  Cross-Section  of  Laundry  Passageway  Section  2 

The  cross  sectional  area  was  calculated  as  a  function  of  the 
coordinates  of  the  points  P0,  P,  and  P2.  The  details  of  this 
calculation  were  discussed  in  reference  [7]. 


4.1.3  Estimation  of  actual  surface  areas 

As  in  the  case  of  the  Laundry  Room,  the  Laundry  Passageway  bulkheads  were  made  of  several 
different  thicknesses  of  steel  plate.  We  again  created  a  new  material  for  the  thermophysical 
database  that  had  a  thickness  equal  to  the  area-weighted  mean  thickness  of  the  actual  materials. 
In  order  to  calculate  this  average,  we  first  estimated  the  areas  for  each  portion  of  the  real 
bulkheads. 

For  Section  1,  these  areas  were  immediately  found  from  the  dimensions 


^dcck  —  AmW  —  4.27  in  - 

Ecjn.  13a 

A,w.i  =  A,„  =  4.50  m: 

Eqn.  13b 

A,lW=6.27  m: 

Eqn.  13c 

Note  that  the  entire  port  side  of  Section  1  is  open  to  Section  2  and,  therefore,  there  is  no  A|loll. 

In  the  case  of  Section  2,  the  areas  of  the  forward  and  aft  bulkheads  were  calculated  from  our 
cross-section  area  expression,  A(w),  using  the  appropriate  values  for  w.  The  slanted  portion  of 
the  port  bulkhead  area  was  treated  as  a  trapezoid  having  two  parallel  edees  the  lenelhs  of  which 
were  the  hypotenuses  of  the  dotted  triangles  shown  in  Figure  9.  The  vertical  part  of  (lie  port 
bulkhead  was  also  a  trapezoid,  but  had  to  be  split  into  upper  and  lower  parts  because  the  hull 
thickness  changed  at  the  1.83  in  (6.00  ft)  level. 

Ihus,  the  surface  areas  for  Section  2  were 


At,^k  =  1 3.66  nv 


Hqn.  14a 
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A0Vhd=  19.19  m2 

Eqn.  14b 

Arwd  =  5.90m2 

Eqn.  14c 

Aafl=  5.26  in2 

Eqn.  14d 

As(bd=  15.60  m2 

Eqn.  14e 

Aport  upper  6.30  m2 

Eqn.  14f 

^pon  lower  =  10.00  IT12 

Eqn.  14g 

Apori  slam  ~  7.84  IT12 

Eqn.  14h 

The  deck  snd  overhead  were  made  of  the  samp  0  qs  \  i  • 

itself  (SHIP3/8).  The  area-weighted  mean  ihh?  f  m'37.5  m  iSte^  as  1,1  the  Laundry  Room 

calculated  to  be  1.05  t^a^Y  PassaSeway  bulkhead  was 

each  part  of  each  surface  The  thermal  nronerties  h  h  &  4  *StS  a,reas  and  Sicknesses  for 
SHIPLRP.  Crmal  pi0perties  database  entry  for  this  material  was  labeled 


Surface 

Section  1 

Section  2 

Adeck 

Area  (m2) 

Thick  (cm) 

Area  (m2) 

Thick  (cm) 

4.27 

0.952 

13.66 

0.952 

Aovhd 

4.27 

0.952 

19.19 

0.952 

Afwd 

4.50 

0.318 

5.90 

0.952 

Aalt 

4.50 

0.952 

5.26 

0.952 

Astbd 

6.27 

1.270 

15.60 

0.3 1 8 

Aport.upper 

- 

— 

6.30 

1 .905 

Aport,  lower 

— 

— 

10.00 

1.588 

Aport  ,slant 

— 

— 

7.84 

1.588 

*  able  4-  Laundry  Passageway  Boundary  Areas  and  Thicknesses 

Section  1  has  no  port  bulkhead  and  the  Section  2  port  bulkhead  is 

Slh:^°7"n;Cpar,S  |;‘'r^“vcn,llco„,p,,',„on,  ink,,;: 
lotul  bulkhead  aica  is  meamnglul, 

4.1.4  Approximation  of  compartment  dimensions 

!:■  1°'  ,f  t;-  *■  ■  *■>  »«•  n* 

Since  there  were  no  hi  n  our  v  ,i  f  ’  '  '".  "'f  lm,er  arc  s""  ">  Ik  tic, conned, 

approximate  valuer  lor  h",|!  '  Vjl'";s  lor  C"I,CI  «*'  ''"'tens,,,,,.  we  had  to  use 
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Note  that  most  of  Section  2  is  the  same  height  as  Section  1  (2.57  m)15and  that  all  of  the 
horizontal  vents  were  in  regions  with  that  height.  Recalling  our  discussion  of  the  importance  of 
compartment  and  soffit  heights  for  horizontal  vent  flow,  we  concluded  that  we  could 
approximate  Section  2  as  if  it  all  had  a  height  of  2.57  m  with  minimal  error. 

Accordingly,  we  have 


H,=  2.57  m 


Eqn.  15 


and  it  follows  that 


D2  =  V2/(H,»  W2)  =  2.17  m  Eqn.  16 

Table  5  summarizes  the  approximate  dimensions  that  were  used  for  both  sections  of  the  Laundry 
Passageway. 


Quantity  (units) 

Psgwy  (Sec  1) 

Psgwy  (Sec  2) 

H  (m) 

2.57 

2.57 

W  (m) 

2.44 

8.51 

D  (m) 

1.75 

2.17 

V  (m3) 

10.97 

47.50  j 

Table  5.  Dimensions  of  Laundry  Passageway  Section 


4.2  Geometrical  Simplifications 

Based  on  the  approximate  compartment  dimensions  discussed  above,  tw'o  alternative  approaches 
to  modeling  the  Laundry  Passageway  u'ere  developed  and  tested.  One  approach  was  to  rearrange 
the  two  sections  of  the  compartment  to  create  a  single  rectangular  parallelepiped  which  could  be 
represented  by  one  set  of  dimensions.  The  alternative  method  involved  treating  each  section  of 
the  Laundry  Passageway  as  a  separate  compartment.  We  refer  to  this  as  the  virtual  compartment 
technique  because  it  resulted  in  the  inclusion  of  an  additional  compartment  in  the  model  -  one 
that  does  not  actually  exist. 

4.2.1  Rearranging  the  compartment 

Our  goal  in  rearranging  the  pieces  of  the  Laundry  Passageway  was  to  create  a  set  of  effective 
dimensions  (Hcfl\  Wcff  and  I\- n)  which  reasonably  approximate  the  compartment. 

From  Equation  7.  the  actual  Laundry  Passageway  volume  was  estimated  to  be 

V  =  58.47  up  Eqn.  17 

Since  the  approximate  height  of  Section  2  was  the  same  as  the  actual  height  of  Section  1,  we 
chose  to  use  that  value  as  the  effective  height 

H,„  =  2.57  m  Eqn.  IS 


The  exception  is  the  narrow  strip,  alone  the  port  edge  of  the  compartment,  where  the  slanted  portion  of  the  hull 
reduced  the  effective  height 


12.  Essentially,  we  have  rotated  Section  2?hv  QnX  h°!f  he  arTangement  shown  in  Figure 
wid,h  and  depth  dimensions  am  swapped  Kf&(whaS?J<:d  k  '°  Section  ,  so  thauhe 

then  the  actual  depth  of  Section  1  plus the  actual  ■ vS  nfw  f  ti?er<5q.ulvaIent  compartment  is 
we  obtained  F  e  aaual  w,dth  of  Section  2.  Using  values  from  Table  5, 


D,  +  W,  =  10.26 


and  the  effective  width  was 


Eqn.  19 


geometry  given  in  Listing  6.  combined  Laundry  Room/Laundry  Passageway 


Ihc  OWCPUon  (°  is  Hull,  by  us  ins:  the  HI/I  ,IU|  CI-VON  i 
positions  ol  .he  floors  of  different  compartments.  bywords,  ,t 


i'  possible  lo  specify  (lie  relative 
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3  SHADWELL/688  Laundry  Passageway  (Rearrangement) 
Sim. time  Print  Hist.  Disp.  Copies 


VERSN  3  SHADWELL/688  3 

#  Sim. time  Print 

TIMES  1250  130 

#  Temp.  Press.  Elev. 

TAMB  285.900  101300.  0.000000 

EAMB  286.300  101300.  0.000000 

#  Crr.pt .  1  Cmpt .  2  Cmpt .  3 

#  Laundry  Psgwy  Exterior 

#Floor  elevation 

HI /F  C.00  0.00 

#X  dimer. . 

DEPTH  1.75  10.26 

#Y  dimer. . 

WIDTH  6.07  2.22 

#Z  dimer.. 

HEIGH  2.57  2.57 

#Materials 

CEILI  SHIP3/8  SHIP3/8 

WALLS  SHIPLR  SHIPLRP 

FLOOR  SHIP3/8  SHIP3/8 

# Laundry -Passageway  door 


3  0 

Elev . 
0.000000 
0.000000 


# 

Cm.pt  #  Cmpt# 

Vent# 

Width 

Soffit 

HVENT 

1  2 

1 

0.66 

1.90 

# 

Crr.pt  #  Cmpt# 

Vent# 

Width@t0 

Width@tl 

CVENT 

1*  2 

1 

1.00 

1.00 

#Passagev.*ay-AMR  door 

'# 

Cmpt#  Cmpt# 

Vent# 

Width 

Soffit 

HVENT 

2  3 

1 

0.66 

2.04 

# 

Cmpt#  Cmpt# 

Vent# 

Width@t0 

Width® tl 

CVENT 

2  3 

1 

1.00 

1.00 

#Passagev.'ay-Torpedo  Rm  door 

# 

Cmpt#  Cmpt  It 

Vent# 

Width 

Soffit 

HVENT 

2  3 

2 

0.66 

2.04 

# 

Cmpt#  Cmpt# 

Vent# 

Widtli@t0 

Width@tl 

CVENT 

2  3 

2 

O 

o 

1.00 

#  Passageway-Wardroom 

hatch 

# 

Cmpt#  Cmpt# 

Area 

Type  ( 1  = 

circular 

WENT 

2  3 

0.25 

2 

# 

y  . 

Z 

EPOS 

C . 91  1.83  0. 

19 

#Fire 

Cmpt 

LFBO 

1 

#  Fi  re 

T\*pe  ( 1  -  uncon 

strained 

;  2  -  con 

strained) 

LFBT 

2 

# 

to 

tl 

FTIME 

1250 

#Mass 

pyrolysis  rate 

FMASS 

0.0253 

0.0229 

# 

Me 1  WU  Kol  Hum  LOL 

He  Inii 

CHEMI 

-84.  100.  10. 

4 . 1 9E^ 

007  285 

#H:C  mass  ratio  (fuel 

compos i t ion i 

HCR 

0.143 

0.14 

3 

#  O :  C  ma  s  s  rati o  ( i ue  3 

compos: 

t  i  on ' 

02 

0 . 0 

0 . 0 

Wind 

0.00 


Wind 
0 . 00 


Listing  6.  Cf  AS  I  Input  Kile  lor  the  Laundry  Passageway  Simulation  using  the  Compartment 

Rearrangement  Method 


#Soot:C02  mass  ratio  (combustion)  ’  - - — - - - 

^  0.06  0.06 

#C0:C02  mass  ratio  (combustion) 

C0  0.056  0.056 

#HCN:fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HCl:fuel  mass  ratio  (pyrolysis) 

HCL  0.0  0.0 

#Toxics : fuel  mass  ratio  (pyrolysis) 

CT  0.0  0.0 

CJET  OFF 

DUMPR  model. HI 

Listing  6  (cont'd).  CFAST  Input  File  for  the  Laundry  Passageway  Simulation  usin«  the 
1 - - - - Compartment  Rearrangement  Method _ 

4.2.2  Use  of  virtual  compartments 

rearrangement  merto^Sau^Mch^inualltomK^ent8!^8'  °f  free?.om  lha"  does  lhe 
definitions.  This  permits  the  actual  rehtinn.hm.  P  ‘  h  S  "S  0wn  dlmens'°ns  and  vent 


S  ec  tdo 


S  ec  tdo  n 


Figure  13.  “Virtual  C 


omparlinenl  Representation  of'  the  Laundry 
Passageway 


!Xilm,3DPT^l":C..UV  !llc  dV"0:!sions  IV<""  Table  5  to  define  , 


taken  together,  represent  the  Laundry  Passage 


wo  compartments  which. 


eway.  In  (his  model,  the  lire  compartment  (Laundry 


Room)  was  compartment  one  so  Sections  1  and  2  of  the  Laundry  Passageway  became 
compartments  two  and  three,  respectively. 

In  this  representation,  the  vertical  vent  and  one  of  the  horizontal  vents  is  located  in  compartment 
three  and  the  second  horizontal  vent  is  in  compartment  two.  All  three  connect  to  the  external 
environment,  which  is  now  compartment  four.  Finally,  a  new  horizontal  vent  was  defined  to 
represent  the  entire  2.44  m  (8.01  ft)  x  2.57  m  (8.43  ft)  bulkhead  between  the  two  sections  The 
equivalent  geometry  of  this  layout  is  shown  in  Figure  13. 

The  area- weighted  bulkhead  thicknesses  were  recalculated,  using  the  values  from  Table  4  to  be 
0.0090  cm  (0.353  in.)  and  0.0110  cm  (0.433  in.)  for  compartments  two  and  three,  respectively 
Two  new  entries,  SHIPLRP1  and  SHIPLRP2,  were  defined  in  the  Thermal.df  file  to  represent 
these  materials.  The  resulting  model  geometry  is  shown  in  Listing  7. 


VERSN 

# 

TIMES 

# 

TAMB 

EAMB 

# 

# 


3  SHADWELL/688  Laundry  Passageway  (Virtual) 


Sim. time 
1250 
Temp . 
285.900 
286.300 
Cmpt.  1 
Laundry 


#Floor  elevation 
HI / F  0.00 

#X  dimer. . 

DEPTH  1.75 
#Y  dimer. . 

WIDTH  6 . 07 
#Z  dimen. 

HEIGH  2 . 57 
#Mater ial s 
CEILI  SHIP3/8 
WALLS  SH I PLR 


Print 

1 

Press . 
101300. 
101300. 
Cmpt .  2 
Psgwy  1 


Hist.  Disp. 
3  0 

Elev. 

0 . 000000 
0.000000 
Cmpt .  3 
Psgwy  2 


Copies 

0 


Exterior 


FLOOR  SH I P3/ 8 
# Laundry- Psgwy  1 

#  Cmpt#  Cm 

HVENT  1  2 

#  Cmpt#  Cm 

CVENT  1  *  2 

#  Psgwy  I -AMR  doo 

#  Cmpt#  Cm 

HVENT  2  4 

#  Cmpt#  Cm; 

CVENT  2  4 


HVENT 

# 

CVENT 

#  Psgwy 

# 

HVENT 

# 

CVENT 


SHIP3/8 

SHIPLRP1 

SHIP3/8 

door 


SHIP3/8 

SHIPLRP2 

SHIP3/8 


Cmpt#  Cmpt# 

Vent# 

Width 

Soffit 

Sill 

Wind 

1  2 

1 

0.66 

1.90 

0.00 

0 . 00 

Cmpt#  Cmpt# 

Vent# 

Width@t0 

Width@t 1 

1  '  2 

1 

1.00 

1 . 00 

r  I -AMR  door 

Cmpt#  Cmpt# 

Vent# 

Width 

Soffit 

Sill 

Wind 

2  4 

1 

0.66 

2.04 

0.23 

0.00 

Cmpt#  Cmpt# 

Vent# 

Width® tO 

Width@t 1 

2  4 

1 

1.00 

1 . 00 

r  1- Psgwy  2  door 

Cmpt#  Cmpt# 

Vent# 

Width 

Soffit 

Sill 

Wind 

2  3 

1 

2.44 

2 . 57 

0.00 

0 . 00 

Cmpt#  Cmpt# 

Vent# 

Width® tO 

Width® t 1 

2  3 

1 

1 . 00 

1  .  00 

■  2 -Torpedo  Rm  door 

Cmpt#  Cmpt# 

Vent# 

Widt.  h 

Sof  fit 

Sill 

Wind 

3  4 

1 

0 . 66 

2.04 

0.23 

0 . 00 

Cmpt#  Cmpt  # 

Vent# 

Wi  dt.hOtO 

Wi dth® t 1 

3  4 

1 

1  .  00 

1  .  00 

Listing  7.  C  1  AS  I  Input  Hie  lot  the  Laundry  Passageway  Simulation  using  the  Virtual 
_  Compartment  Method 
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#Passageway-Wardroom  hatch  - - - 

WENT  C3Pt‘  TT  U  =  Circul«''  2 

#  x  Y  z 

FPOS  0.91  1.83  0.19 

#Fire  Cmpt 
LFBO  1 

=  unconstrained;  2  =  constrained) 


square) 


He 

4 . 19E+007 


Init  T 
285.9 


T  Rad.  tract. 
330.  0.30 


#  to  tl 

FTIME  1250. 

#Mass  pyrolysis  rate 

FMASS  0.0253  0.0229 

#  Mol  Wt  Rel  Hum  LOL  He  Tn(t.  „ 

CliEMI  184.  100.  10.  4 . 19E+007  285  9  T  ***'  fraCt‘ 

#H:C  mass  ratio  (fuel  composition)  85’9  °'30 

HCR  0.143  0.143 

#0:C  mass  ratio  (fuel  composition) 

02  0.0  0.0 

#Soot:C02  mass  ratio  (combustion) 

00  0.06  0.06 

#C0:C02  mass  ratio  (combustion) 

C0  0.056  0.056 

#HCN : fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HCl:fuel  mass  ratio  (pyrolysis) 

HCL  0.0  0.0 

KToxics : fuel  mass  ratio  (pyrolysis) 

CT  O.o  0.0 

CJET  OFF 
DUMPR  model. HI 

LiS,i"S  7  <C0",'d)'  CFAST  Laundry  Passageway  Simula, ion  using  ,„e  Vim,,,! 

1 - - - - - Compartment  Method _ ^  lluial 

4.3  Problems  of  Vertical  Transport 

-ss  transport  through  vertical 

aspect  was  first  encountered  when  the  ve^i  i  t  J„ ,  Y ,adja“nt  compartments.  The  former 
the  model.  C  Ve,l"‘jl  ve,u  1,1  ,he  Sundry  Passageway  was  included  in 

™yvi™C^&STa,'^  boundary  is  exposed  ,hc  ex, on, a, 

compartment,  is  used  for  conduction  calcul  uiorw  li  A’  'J,  ,CI  1  ldn  lhe  tcmPe|ature  of  another 
external  environment  i„sS of in  „  ‘he seZ)l  ‘'l?*  ener&  is  'ransferred  to  the 

model  this  behavior  was  coVrncl  -  the  far  of  ,U  T'l '"T .  KT  lhe  Lai,,,dl>  Passageway 

the  Wardroom  was  added  comparlmenl-lo  eomn-in  °  >tn  cI  Katl  the  exterior.  However,  when 
the  default  behavior  was  no  longer  appropriate.  vcrtical  conduction  became  a  factor  and 

Mh  of  ,l,ese  vertical  ,ransp„r, 
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4.3.1 


Vertical  mass  transport 


The  Laundry  Passageway  vertical  vent  was  rectangular,  with  dimensions  of  0.96  m  (3.15  ft)  by 
0.81  m  (2.66  ft).  Because  CFAST  limits  VVENTs  to  circles  or  squares,  we  had  to  approximate 
this  as  a  square,  which  was  reasonably  close  to  the  actual  shape.  The  vent  area  was  defined  to  be 
equal  to  that  of  the  actual  hatch,  0.78  m2  (8.40  ft2). 

The  hatch  described  above  was  included  in  both  the  "rearranged  geometry"  and  the  "virtual 
compartment"  versions  of  the  Laundry  Passageway  model.  In  both  cases,  we  found  that  the 
model  failed  to  run  to  the  end  of  the  1250  second  simulation.  When  the  "rearranged  geometry" 
approach  was  used,  the  model  ran  at  normal  speeds  for  the  first  148  seconds  of  the  simulation  (at 
a  rate  of  about  four  simulated  seconds  per  second  of  computer  time17)  but  then  the  execution 
speed  dropped  four  orders  of  magnitude  (to  approximately  3  simulated  milliseconds  per  second). 
With  the  virtual  compartment  model,  CFAST  stalled  after  only  78  seconds  of  simulation.  At 
those  rates,  the  model  is  effectively  unusable  because,  on  typical  desktop  computers,  it  would 
take  days  to  complete  each  simulation. 

The  stalling  problem  is  a  consequence  of  the  way  in  which  CFAST  calculates  the  size  of  the  time 
step  to  be  used  by  the  numerical  solver.  The  step  size  is  adjusted,  upward  or  downward,  based  on 
the  results  of  previous  steps.  The  step  size  is  increased  when  convergence  is  reached  very 
quickly  and  decreased  if  convergence  fails.  Typically,  the  step  size  is  on  the  order  of 
milliseconds  (it  can  be  as  much  as  several  seconds  when  conditions  are  especially  favorable)  but, 
after  repeated  convergence  failures,  the  step  size  can  be  reduced  to  a  fraction  of  a  microsecond. 
Under  normal  circumstances,  simulating  one  second  requires  only  hundreds  or  thousands  of 
iterations  but  this  can  grow  to  tens  of  millions  of  iterations  under  adverse  conditions. 


When  faced  with  stalling  problems,  the  first  step  is  to  determine  which  keywords  are  involved1*. 
The  easiest  way  to  accomplish  this  is  to  "comment  out”  the  suspect  lines,  one  at  a  time,  by 
inserting  a  pound  sign  (#)  as  the  first  character  in  the  line.  If  the  model  no  longer  stalls,  then 
some  parameter  associated  with  that  keyword  was  a  contributing  factor. 


Assuming  that  WENT  has  been  so  identified,  there  are  several  possible  methods  for  resolving 
the  problem.  The  first  is  to  try  different  vent  areas  or  shapes  until  a  combination  is  found  which 
permits  the  model  to  run  to  completion.  It  is  then  the  responsibility  of  the  user  to  determine 
whether  the  resulting  parameters  adequately  represent  the  actual  situation.  If  they  do  not,  then  the 
only  iecour.se  is  to  identify  additional  interacting  keywords  or  parameters  and  make  adjustments 
to  them. 


Using  the  compartment  rearrangement  method,  we  found  that  the  Laundry  Passageway  model 
ran  for  the  full  simulation  period  if  the  area  of  the  WENT  was  less  than  about  0.25  m2  (2.69  ft2) 
or  greater  than  approximately  7.5  in2  (66.2  ft2).  With  the  virtual  compartment  technique,  serious 
stalling  problems  were  noted  for  all  non-zero  vent  areas.  These  results  are  illustrated  in  Figures 
14  and  15.  which  show  the  completion  fraction,  as  a  function  of  vent  area,  for  the  two  different 


Typical  execution  rates  on  current  microcomputers  are  on  the  order  of  It)  simulated  seconds  per  second.  However, 
the  solver  is  slower  at  the  beginning  of  the  simulation  because  it  is  making  essentially  random  guesses.  Once  it  has 
acquired  a  history  of  prior  solutions,  the  solver  predictions  gel  better  and  it  speeds  up.  Transient  inputs  (a  step 
function  change  in  burning  rate,  (or  example)  normally  cause  a  temporary  slowdown  until  a  new  history  can  be 
developed  Also.  ('FAST  requires  several  minutes  to  load  and  initialize.  Since  this  disproportionately  affects  short 
simulations,  we  did  not  begin  timing  until  after  the  startup  delay. 

Runtime  pioblcms  with  (  FAS  I  aie  usually  due  to  interactions  among  parameters  associated  with  multiple* 
keywords  1  here! ore,  it  is  not  generally  correct  to  say  that  a  specific  parameter  or  keyword  is  the  cause  of  the 

difficulties 
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4 


It  became  apparent  that,  if  we  used  the  virtual  compartment  method,  we  would  have  to  eliminate 
the  Laundry  Passageway-Wardroom  hatch  from  the  model.  Since  it  was  obvious  that  the  model 
predictions  for  the  Wardroom  would  be  grossly  incorrect  without  that  hatch,  we  chose  not  to 
pursue  the  virtual  compartment  approach  any  further  for  this  work.  We  must  note,  however  that 
the  tact  that  virtual  compartments  were  inappropriate  in  this  specific  case  does  not  indicate  that 
they  could  not  be  useful  in  other  cases. 

We  have  found  that  stalling  problems  sometimes  vanish  if  the  dimensions  of  one  of  the 
comPfrtments  1S  changed.  For  example,  our  difficulties  with  the  Laundry  Passageway 
V V ENT  were  eliminated  when  the  Wardroom  was  added.  The  effect  of  this  addition  was  to 
change  the  destination  compartment  from  the  essentially  infinite  volume  of  the  exterior  to  a  finite 

V  Ui  Li  1  ilC- 

Adding  the  Wardroom  might  also  have  made  it  possible  to  use  virtual  compartments 
Unfortunately,  there  are  so  many  possible  approximations  for  complex  models  that  combinatorial' 
explosion  becomes  a  problem  —  it  is  impossible  to  exhaustively  investigate  all  combinations  and 
the  user  must  prune  entire  branches  from  the  possibility  tree.  In  our" case,  we  elected  not  to 
pursue  the  virtual  compartment  approach  and,  instead,  used  compartment  rearrangement.  The 

vent  area  was  set  to  0.25  m2  (2.69  ft2)  because  that  was  the  closest  we  could  get  to  the  actual 
hatch  area. 


4.3.2  Vertical  heat  conduction 

By  default,  CFAST  does  not  know  that  the  overhead  of  one  compartment  is  the  deck  of  another. 
i  nerefore  heat  conduction  between  two  vertically  adjacent  compartments  is  normally  not 
calculated.  However,  the  CFCON  keyword  can  be  used  to  inform  CFAST  that  a  particular  pair  of 
compartments  is  vertically  adjacent.  When  used,  this  keyword  enables  vertical  heat  transfer  from 
the  lower  compaitment  to  the  upper19.  Modeling  of  the  Wardroom  provided  an  example  of  the 
application  of  vertical  heat  conduction.  F 

Since  the  Wardroom  overlies  both  the  Laundry  Room  and  the  Laundrv  Passageway  the  obvious 
approach  was  to  connect  both  lower-deck  compartments  to  the  Wardroom.  Unfortunately  this 
did  not  work  —  if  two  CFCON  inputs  to  the  same  compartment  are  specified,  the  numerical 
solver  fails  to  initialize  properly  and  causes  CFAST  to  immediately  quit2".  Consequently  in 
cases  like  this,  the  user  must  determine  which  one  of  the  possible'  connections  is  the  most 
important  and  which  can  salely  be  neglected. 

We  investigated  both  possibilities:  conduction  between  Laundry  and  Wardroom  and  between 
Laundry  Passageway  and  Wardroom.  There  was  an  interaction  between  the  Wardroom-exterior 
Sm^E|NT-  ,a,ld  thC  La  11 " d ry - w a rd roo in  or  Laundry  Passageway-Wardroom  conduction 
(CFCON).  In  either  case,  when  the  correct  area  10.78  m2  (8.40  ft2)]  of  the  exterior  vertical  vent 
was  specified,  there  were  stalling  problems  similar  to  those  previously  discussed  Without 
vertical  conduction,  the  model  ran  correctly  with  the  actual  vent  area.  By  trial  and  error  it  was 
discovered  that  the  model  would  run  for  either  conduction  configuration  if  the  exterior  vent  size 
was  reduced  to  0.26  m2  (2.80  ft2)  or  below.  '  ‘ 

We  chose  to  use  conduction  between  the  Laundry  Room  and  the  Wardroom  for  two  reasons:  1) 
conduction  Irom  the  Laundry  was  greater  than  conduction  from  the  Laundry  Passageway 

jq 

(ICON  changes  Hie  behavior  of  vertical  conduction  (through  overheads/decks)  to  permit  enercv  to  he  added  to 
die  uppei  compartment.  Ii  does  not  ailed  horizontal  heat  conduction  (through  bulkheads) 

j  h  ,S  :'ccci;,ablc  ,or  ono  compartment  to  have  both  a  conductive  input  (through  the  deck)  and  an  output  (throu-h 
the  overhead  I.  1  r 
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froT,h=‘^ 


1250 
Temp . 
285.900 
286.300 
Cmpt .  1 
Laundry 
#Floor  elevation 
HI/F  0.00 
#X  dimen . 

DEPTH  1.75 
#Y  dimen. 

WIDTH  6.07 
#Z  dimen. 

HEIGH  2.57 
#Materials 
CEILI  SHIP3/8 
WALLS  SHIPLR 
FLOOR  SHIP3/8 


3  SHADWELL/ 688  Laundry  -  Wardroom 
Sim. time  Print  Hist.  Disp.  Copies 

3  0 


1 

Press . 
101300. 
101300. 
Cmpt.  2 
Psgwy 

0.00 

10.26 

2.22 

2.57 

SHIP3/8 

SHIPLRP 

SHIP3/8 


Elev. 
0.000000 
0.000000 
Cmpt .  3 
Wardroom 

2 . 57 


0 


Cmpt .  4 
Exterior 


#Laundry— Passageway  door 

4  _ .  „  _ 


Vent# 


#  Cmpt# 

HVENT  1 

#  Cmpt#  Cmpt#  Vent# 

CVENT  1  2  1 

# Pas sageway- AMR  door 

#  Cmpt#  Cmpt#  Vent# 

HVENT  2  4  x 

#  Cmpt#  Cmpt#  Vent# 

CVENT  2  4  1 

# Passageway-Torpedo  Rm  door 

#  C^ot#  Cmpt#  vent# 

HVENT  242 

#  Cmpt#  Cmpt a  Vent# 

CVENT  242 

# Passageway-Wardroom  hatch 

#  Cmpt#  Cmpt#  Area 

WENT  2  3  0.78 

# Ward room -Crew  Mess  door 

Cmpt#  Cmpt#  Vent# 

HVENT  3  4  2 

#  Cmpt#  Cmpt#  Vent# 

CVENT  34X 

# Wardroom -Crew  Living  door 

#  Cmpt#  Cmpt#  Vent# 

HVENT  342 

#  Cmpt#  Cmpt#  Vent# 

CVENT  349 


4 . 00 

8.51 

2.59 

SHIP7/8 

SHIPWR 

SHIP3/8 

Width 
0 . 66 

Width@t0 
1 . 00 

Width 
0 . 66 

Width@t0 
1 . 00 

Width 
0 . 66 

Width© to 
1.00 


# Ward room - Na v .  Equip. 

#  Cmpt#  Cmpt# 

WENT  3  4 

#  Cmpt#  Cmpt# 


Width 

0.66 

WidthOt 0 
1.00 

Width 
0 . 66 

Width© t0 
1 . 00 


Room  hatch 
Area  Type 
0.26  2 


Soffit 

2.04 

Width@t 1 
1 .00 

Soffit 
2 . 04 

Width@tl 
1  .  00 


Soffit 
1 .90 

Width© tl 
1.00 

Soffit 

2.04 

Width© tl 
1.00 

Soffit 

2.04 

Width@t 1 
1.00 


Sill 

0.00 


Sill 

0.23 


Wind 
0 . 00 


(1  -  circular;  2  = 


square) 


Sill 

0.23 


)  .23 


Wind 
0 . 00 


tl 


circular;  2 


square ) 


-Uslino  S.  Cl-'AST  lupin  l-ile  for  the  Wardroom  Simula! 


ion 
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CFCON  1  3 

#  X  Y  Z 

FPOS  0.91  1.83  0.19 

#Fire  Cmpt 
LFBO  1 

#Fire  Type  (1  =  unconstrained;  2  =  constrained) 
LFBT  2 

#  tO  tl 

FTIME  1250. 

#Mass  pyrolysis  rate 
FMASS  0.0253  0.0229 

#  Mol  Wt  Rel  Hum  LOL  He 

CHEMI  184.  100.  10.  4.19E+007 

#H:C  mass  ratio  (fuel  composition) 

HCR  0.143  0.143 

#0:C  mass  ratio  (fuel  composition) 

02  0.0  0.0 

#Soot:C02  mass  ratio  (combustion) 

OD  0.06  0.06 

#C0:C02  mass  ratio  (combustion) 

CO  0.056  0.056 

#HCN:fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HCl:fuel  mass  ratio  (pyrolysis) 

HCL  0.0  0.0 

#Toxics:fuel  mass  ratio  (pyrolysis) 

CT  0.0  0.0 

CJET  OFF 
DUMPR  model . HI 


Init  T 
285.9 


Ign.  T  Rad.  f ract . 
330.  0.30 


_ Listing  8  (cont'd).  CFAST  Input  File  for  the  Wardroom  Simulation _ 

4.4  Problems  of  Ambiguous  Thermophysical  Properties 

We  have  seen,  in  both  the  Laundry  Room  and  (he  Laundry  Passageway,  that  it  is  necessary  to 
create  fictitious  materials  in  cases  where  the  actual  bulkheads  are  composed  of  patches  of 
different  materials21.  In  the  previous  cases,  this  problem  was  relatively  simple  because  the 
patches  had  the  same  composition  and  differed  only  in  thickness.  The  Navigation  Equipment 
Room  (Figure  16)  was  much  more  complex  because  the  bulkheads  differed  in  composition  as 
well  as  in  thickness. 

All  of  the  original  bulkheads  in  the  Navigation  Equipment  Room  were  composed  of  various 
thicknesses  of  steel.  However,  the  forward  bulkhead  (on  the  right  in  Figure  16)  was  a  temporary 
partition  made  of  plywood.  The  thermal  conductivities,  heal  capacities  and  densities  of  steel  and 
plywood  arc  vastly  different  and  it  was  not  immediately  obvious  what  weighting  factors  should 
be  used  to  calculate  meaningful  average  values. 

For  consistency  with  our  prior  work,  we  used  area-weighting  to  calculate  the  mean  thickness  and 
this  thickness  was  used  for  all  subsequent  calculations.  For  the  conductivity,  we  required  that  the 
effective  heat  flux.  Q,  be  equal  to  the  sum  of  the  actual  heat  fluxes.  Q,.  through  the  various  pieces 
of  the  bulkhead 


Q  =  £  Q, 


Eqn.  21 


This  rule  also  applies  to  deck  and  overhead  boundaries,  although  there  were  no  such  examples  in  this  ease  study. 
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or,  in  terms  of  the  material  properties, 


kAaT  / 1  =  2)  KjAj  ATj  / 1| 


Eqn.  22 


the  thickness.  The  s ubscripts’ refer  t^h^eSionronhe^Ikh  temperature  difference  across 
the  effective  values  for  our  approximated  compartment  unsubscriP'«i  values  refer  ,o 


k=Zk,(A,/A)(i/(,) 


Iiqn.  23 


thickness  assumed 

“ . . . -• .  and  i:;:;'. 


tliai  the  total  mass  of  the 
bulkheads.  Applying  this 
ives  us 
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p  =  2  Pi(Aj  /  A)  (tj  /  t) 


Eqn.  24 


Similarly,  the  requirement  for  equal  heat  content  of  the  bulkheads  leads  to 


C  =  ECi(pi/p)(Ai/A)  (t,/t)  Eqn.  25 

In  order  to  apply  these  equations,  we  first  calculated  the  areas  of  the  seven  parts  of  the 
Navigation  Equipment  Room  bulkhead  (indicated  by  the  circled  numbers  in  Figure  17)  using  the 
compartment  dimensions  from  Figure  16.  The  area- weighted  thickness  was  then  calculated  and, 
with  this  information,  we  were  able  to  determine  all  of  the  weighting  factors  for  Equations  23  - 
25.  Standard  values  for  the  thermophysical  properties  of  steel  and  plywood  were  taken  from  the 
default  Thermal.df  file  (the  relevant  portion  of  which  is  shown  in  Table  6).  The  calculated  mean 
property  values  for  the  fictitious  material  (SHIPNER)  are  given  in  Table  7. 


Figure  17.  Division  of  the  Navigation  Equipment  Room  Deck  into  Parts 

The  deck  area  was  calculated  as  the  sum  of  the  areas  of  the  three  parts 
illustrated.  In  order  to  estimate  mean  values  for  the  thermophysical 
properties,  the  bulkhead  was  divided  into  seven  segments,  as  indicated  by 
the  circled  numbers.  The  area  of  each  segment  was  the  length  of  the  line 
segment  multiplied  by  the  compartment  height. _ 


Material 

Conductivity 

Heat  Capacity 

Density 

(W/m  K) 

(J/kg  K) 

(kg/m3) 

Steel 

48 

559 

7854 

Plywood 

0.12 

1215 

545 

Table  6.  Standard  Thermal.df  Entries  for  Steel  and  Plvwood 
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1 

2 

3 

4* 

5* 

6* 

7 

Totals 


5.23 
3.81 
3.29 
4.17 
3.78 
1.19 

3.24 
24.71 


0.0095 

0.0095 

0.0095 

0.0127 

0.0127 

0.0127 

0.0095 


.2117 

.1542 

1331 

1688 

1530 

0482 

1311 


Table  7.  Calculation  of  Thermal  Properties  for  SHIPNER 

*£,e 


'*  °™2  11.34  r55o - 777 - 

f  0K)«  8.26  ,085  no 

2*  0.8962  7.13  937  ,, 

4  1-1981  0.0169  no  29‘ 

5  11981  001 63  100  o6 

6  1-1981  0.0048  3i  “ 

7*  0.8962  7J12  993  lm 

I— I  I  33.79  I  4676  | 

* - — ble  7  (Cont  d).  Calculation  of  Thermal  Properties  for  SHIPNER 

4.5  Results  of  the  Geometry  Specification  Development 

10  data  from  the  SHADWELL/688  ' 'tesl1 1‘ G?1  'l n  t h  ' D  'wT'  °l  ^  ^  Spcal k'a,lon- 
portion  of  the  model  that  was  directly  affected  hv  ii  ,  f;  >  ,  rela,lvel>  easY  «»  separate  the 

Hie  portions  that  were  not  (everythimi  else)  H  )w[w^  (lhe  L;U"Uliy  Roo"»  horn 

interactions  among  the  various  parts  of  the  n  J,el  “£ ‘hc 

i„  . . 

!Kl‘k'11 ,ho  ,mHiei  tiuk-  a"h™8" « i-c 

I'., .h,iiSc.l  if  , .1 'y ' h '' h.,  M I  !m ,h.'i h  'l'.'! ' " . .  '  "  'V  sl»'"  "■  a.I.I.ii..,,  „f  ,|,c 


geometry,  vertical  transport  and  thermophysical  properties  as  if  they  were  independent,  in  fact 
they  are  not  independent.  Accordingly,  in  this  section,  we  compare" the  results  of  modeling  the 
entire  domain  (Listing  9),  rather  than  subsections  of  the  domain,  with  the  test  data.  Even  then 
the  predictions  can  only  be  considered  to  be  estimates  because  our  model  domain  did  not  include 
the  entire  SHADWELL/688  test  area.  As  is  likely  to  be  true  for  most  real  world  applications,  the 
model  domain  was  chosen  to  include  the  factors  of  primary  interest  (in  this  case,  the  effects  of 
Laundry  Room  fires  on  Control  Room  habitability)  while  neglecting  other  factors  (such  as  the 
effects  on  the  AMR,  Torpedo  Room  and  other  compartments). 


VERSN 

3  SHADWELL/688 

Laundry  -  Sail. 

# 

Sim. time 

Print 

Hist.  Disp.  Copies 

TIMES 

1250 

1 

3  0  0 

# 

Temp . 

Press . 

Elev. 

TAMB 

285.900 

101300. 

0.000000 

EAMB 

286.300 

101300 . 

0.000000 

#  Crept .  1  Cmpt .  2 

#  Laundry  Psgwy 
#Floor  elevation 

Cmpt.  3 
Wardrm 

Cmpt .  4 
NER 

Cmpt.  5 
CR 

Cmpt .  6 
Sail_i 

Cmpt .  7 
Sail_2 

HI  / F  C.00 

#X  dimer. . 

0 . 00 

2.57 

5.16 

5.16 

7.75 

10.18 

DEPTH  1.75 
#Y  dimer. . 

10.26 

4 . 00 

2 . 92 

2.90 

1.21 

0.91 

WIDTH  c . 07 
#Z  dimer.. 

2.22 

8.51 

1.90 

6.30 

1.42 

0 . 91 

HEIGH  2 . 57 

2 . 57 

2.59 

2.59 

2 . 59 

2.43 

1 . 02 

^Material s 

CEILI  SHIP3/8 
WALLS  SH I PLR 

FLOOR  SH I P3/ 8 

SHIP3/8 

SHIPLRP 

SHIP3/8 

SHIP7/8 

SHIPWR 

SHIP3/8 

SHIP3/8 

SHIPNER 

SHIP7/8 

SHIP3/8 

SHIPCR 

SHIP7/8 

SHIP3  /  8 
SHIPS  /  8 
SHIPS  -'8 

SHIP3/8 

SHIP3/8 

SHIP3/8 

# Laundry- passageway  door 


# 

C-?:#  Cmpt d 

Vent# 

Width 

Soffit 

Sill 

‘Wind 

HVENT 

:  2 

1 

0 . 66 

1.90 

0 . 00 

o 

o 

o 

# 

c-.pt  if  Cmpt# 

Vent# 

Width® tO 

Width@t 1 

CVENT 

i  2 

1 

1.00 

1.00 

#  Pa  ssaa  ev.*a  y  -  AMR  door 

# 

C.T.p  L  #  Cmpt# 

Vent# 

Width 

Sof  f i t 

Sill 

Wind 

HVENT 

2  8 

1 

0.66 

2.04 

0.23 

0 . 00 

# 

Crr.pt:  #  Cmpt# 

Vent# 

Width@t0 

Width® t 1 

CVENT 

2  8 

1 

1 .00 

1  .00 

#  Pa  ssagevray-  Torpedo 

Rm  door 

# 

Cmpt#  Cmpt.# 

Vent  # 

W  i  d  t  h 

Soffit 

Sill 

Wind 

HVENT 

2  8 

2 

0.66 

2 . 01 

0.23 

0.00 

# 

Crept#  Cmpt  # 

Vent# 

Wi dth® tO 

Wi  dt  h@ 1. 1 

CVENT 

2  8 

2 

1  .00 

1  .  00 

#  Passageway  -  Wa  rdroom 

hatch 

# 

Crr.pt  #  Cmpt  # 

Area 

Type  ( 1  = 

ci rcular ; 

2  =  sc 

rare) 

WENT 

2  3 

0 . 78 

2 

#Wardr 

' o  r  r.  -  C  i  ew  Mess 

door 

tt 

C:r.::t#  C'mpt  # 

Vent# 

Width 

Sof  f i t 

Si  I  1 

■-*  -  -s 

HVENT 

3  8 

1 

0.66 

2 . 04 

0 . 2  3 

o .  o : 

Listing  9.  Hie  lor  the  Hnal  SH  AD\YULL/6tXS  Simulation 
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#  Cmpt#  Cmpt#  Vent#  WidthOtO 

WENT  3  8  1  i.oo 

#Wardroom-Crew  Living  door 

#  Cmpt#  Cmpt#  Vent#  Width 

WENT  3  8  2  0.66 

#  Cmpt#  Cmpt#  Vent#  Width@tO 

CVENT  3  8  2  1.00 

#Wardroom-Nav .  Equip.  Rm.  hatch 

#  Cmpt#  Cmpt#  Area  Type  (1  = 

WENT  3  4  0.78  2 

#Laundry-Wardroom  heat  conduction 

#  Cmpt#  Cmpt# 

CFCON  1  3 


3  Width@tl 
1.00 

Soffit  Sill  wind 

2-04  0.23  0.00 

)  WidthQtl 
1.00 

=  circular;  2  =  square) 


#Nav.  Equip.  Rm-Fan  Rm.  door 
#  Cmpt#  Cmpt#  Vent#  Width 

WENT  4  8  1  0.66 


#  Cmpt#  Cmpt#  Vent#  Width@t0 

CVENT  _4  8  1  i.oo 

#Nav.  -quip.  Rm. -Control  Rm.  door 

#  Cmpt#  Cmpt#  Vent#  Width 

WENT  4  5  1  o.50 

#  Cmpt#  Cmpt#  Vent#  Width@t0 

CVENT  4  6  i  i  ™ 


Soffit  Sill 

2.04  0.23 

Width@tl 
1.00 

Soffit  Sill 

1-92  0.00 

Width@tl 
1.00 


#Contrc-_  Rm.  Combat  Systems  door 

#  Cmpt#  Cmpt#  Vent#  Width 

WENT  5  8  i  o.66 

#  Cmpt#  Cmpt#  Vent#  Width@t0 

WENT  5  8  1  i.oo 

#Contrc 1  Rm.-Sail_l  hatch 

#  Cmpt#  Cmpt#  Area  Type  (1  - 

WENT  5  6  0.33  i 

# Wardroom- Control  Rm.  heat  conduction 

#  Cmpt#  Cmpt# 

CFCON  3  5 


Soffit 
2 . 04 

Width@  tl 
1 . 00 


Wind 

0.00 


Wind 

0.00 


Sill  Wind 

0.23  0.00 


circular;  2  =  square) 


^ Sai 1__ - Sa i 1_2  opening 
^  Cup  tfl  Cmpt#  Area  Type  (1 

WENT  c  7  0.82  1 

ilContrc*  Rm.-Sail_l  heat  conduction 
#  Crr.p  t#  Cmpt  # 

CFCON  5  6 


=' circular;  2  =  square) 


#Sa exterior  opening 

Cu.pt  #  Cmpt#  Area  Type  (1 
WENT  ^  8  0 .82  ] 

#Sai 1_. -Sai heal  conduction 
#  Cu.pt  #  Cmpt  H 

CFCON  6  7 

FPOS  :.91  1.83  0.19 


circular;  2  =  scudre1 


#Fir<?  Type  ( 
i  j*’bt  ;> 


um-onsi  rainod;  2  -  ron.si  r<unod) 


—  mS  - tit  tlUlr  IjilflL  SHADWI-U76SS  Simulation 


4X 


#  to  tl  ~  — 

FTIME  1250. 

#Mass  pyrolysis  rate 

FMASS  0.0253  0.0229 

#  Mol  Wt  Rel  Hum  LOL  He  Init  T  Icr. .  T  Rad.  fract. 

CHEMI  184.  100.  10.  4.19E+007  285.9  330.  0.30 

#H:C  mass  ratio  (fuel  composition) 

HCR  0.143  0.143 

#0:C  mass  ratio  (fuel  composition) 

02  0.0  0.0 

#Soot:C02  mass  ratio  (combustion) 

OD  0.06  0.06 

#CO : C02  mass  ratio  (combustion) 

CO  0.056  0.056 

#HCN : fuel  mass  ratio  (pyrolysis) 

HCN  0.0  0.0 

#HC1 : fuel  mass  ratio  (pyrolysis) 

HCL  0.0  0.0 

#Toxics : f uel  mass  ratio  (pyrolysis) 

CT  0.0  0.0 

CJET  OFF 
DUMPR  model .HI 

_ Listing  9  (cont'd).  File  for  the  Final  SHADWELL/688  Simulation _ 

The  primary  purpose  of  this  section  is  to  provide  some  insights  regarding  the  accuracy  of  CFAST 
predictions.  It  is  important  to  note  that  good  agreement  (or  lack  thereof! in  our  test  case  does  not 
guarantee  that  all  simulations  will  be  equally  accurate  (or  inaccurate).  Data  from  the  various 
compartments  were  processed  as  previously  described  and  the  layer  assignments  for  the 
thermocouples  were  again  based  on  the  interface  heights  predicted  by  CFAST.  Since  it  was  not 
practical  to  apportion  sensors  on  a  second-by-second  basis,  we  used  the  mean  of  the  interface 
heights  at  60  and  1250  seconds.  These  times  were  arbitrary,  but  approximately  represented  the 
period  over  which  the  interface  height  was  nearly  constant.  Table  S  shows  the  number  of 
thermocouples  for  both  layers  of  each  compartment  and  the  interface  heights  for  all 
compartments-3  are  shown  in  Figure  18. 


Cinpt. 

Interface  Hi 
(m) 

Upper  TCs 

Lower  TCs 

Laundry 

0.56 

4 

0 

Passageway 

1.54 

2 

4 

Wardroom 

1.64 

9 

20 

Nav.  Equip. 

0.86 

4 

-) 

Control  Rm. 

1.33 

6 

i 

Table  8.  Interface  Heights  and  Apportionment  of  Thermocouples 


I  he  number  of  thermocouples  in  each  layer  ol  each  compartment 
was  determined  based  on  the  interface  heights  predicted  bv 
CFAST. 
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Figure  18.  Compartment  Interface  Heights 

These  CFAST-predicted  interface  heights  were  used  to  determine  which 
|thumocouples  were  in  the  upper  or  lower  layers. _ 
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-  CFAST  -  Expt'l 


Hie  error  bars  represent  one  standard  deviation  in  the  test  data.  The 
Assignment  of  thermocouples  to  the  upper  or  lower  layer  was 
based  on  (he  interface  height  estimated  by  CFAST. _ 
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Figure  20A.  Predicted  versus  Measured  Laundry  Passagew  ay 
Upper  Layer  Temperatures 

The  error  bars  represent  one  standard  deviation  in  the  test  data.  The 
Assignment  of  thermocouples  to  the  upper  or  lower  layer  w'as 
based  on  the  interface  height  estimated  by  CFAST. 
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Hguie  20B.  Predicted  versus  Measured  Laundry  Passageway 
Lower  Layer  Temperatures 


I  he  ei i oi  bats  repiesent  one  standard  deviation  in  the  test  data  Tiie 
Assignment  of  thermocouples  to  the  upper  or  lower  layer  was 
based  on  the  interlace  height  estimated  by  CFAST. 
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Figure  21  A.  Predicted  versus  Measured  Wardroom  Upper  Layer 

Temperalures 

The  error  bars  represent  one  standard  deviation  in  the  test  data.  The 
Assignment  of  thermocouples  to  the  upper  or  lower  layer  was 
based  on  the  interface  height  estimated  by  CFAST. 
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Figure  2 IB.  Predicted  versus  Measured  Wardroom  Lower  Layer 

Temperatures 


The  error  bars  represent  one  standard  deviation  in  the  test  H-ih  ti,  , 
Assignment  of  thermocouples  to  the  upper  or  Ewer 
based  on  the  interlace  height  estimated  by  CFAST. 
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Figure  23B.  Predicted  versus  Measured  Control  Room  Lower- 
Layer  Temperatures 


The  error  bars  represent  one  standard  deviation  in  the 
Assignment  of  thermocouples  to  the  upper  or  low 
based  on  the  interlace  height  estimated  by  CFAST. 


test  data.  The 
er  layer  was 


There  is  reasonable  agreement  between  model  and  test  for  the  upper  layer  of  the  Navigation 
Equipment  Room  but,  given  the  major  discrepancy  in  the  upper  layer  temperatures  for  the 
Wardroom,  this  result  may  be  coincidental.  The  observation  that  the  predicted  lower  layer 
Navigation  Equipment  Room  temperatures  are  only  slightly  higher  than  the  pre-fire  temperatures 
supports  this  conclusion.  Clearly,  a  column  of  superheated  gas  (from  the  upper  layer  of  the 
Wardroom)  could  not  rise  through  the  lower  layer  of  the  Navigation  Equipment  Room  without 
heating  the  latter.  The  situation  in  the  Control  Room  is  similar  —  there  is  no  significant  increase 
in  the  predicted  lower  layer  temperature,  resulting  in  very  poor  agreement  with  experiment  but 
the  agreement  with  the  data  for  the  upper  layer  is  surprisingly  good. 

5.0  EXAMPLES  OF  ALTERNATIVE  INPUTS 

In  the  two  previous  sections,  we  outlined  an  approach  to  building  the  "best"  model  of  a  specific 
test  case  (/.£.,  one  based  on  the  most  accurate  available  inputs).  In  the  process,  we  have  seen 
several  areas  in  which  alternative  sets  of  inputs  were  possible  and  in  which  the  decision  as  to 
which  set  should  be  used  was  somewhat  arbitrary.  In  this  section,  we  will  discuss  the  results 
obtained  when  some  of  those  alternatives  are  used.  Our  purpose  is  to  illustrate  the  general 
magnitudes  of  the  effects  that  might  be  expected. 

5.1  Fire  Specification  Alternatives 

In  the  course  of  developing  the  fire  specification,  we  were  able  to  estimate  values  for  the  fuel 
parameters  based  on  a  priori  knowledge  regarding  the  type  of  fuel.  However,  we  were  not  able 
to  do  the  same  for  the  pyrolysis  and  combustion  parameters  because  they  are  dependent  on  the 
dynamics  of  the  burning  process.  Typically,  the  user  would  make  "best  guess"  estimates  of  the 
nominal  values  for  these  parameters  and  then  bracket  those  with  high  and  low  extremes  In  this 
section,  we  provide  examples  of  the  effects  of  some  of  the  kev  pyrolysis  and  combustion 
parameters. 


5.1.1  Heat  release  rate  parameters 

Equation  1  showed  the  relationship  among  AHC,  FMASS  and  FQDOT,  the  three  parameters 
relevant  to  heat  release  rale.  In  our  case  study,  we  specified  the  first  two  of  these  and  CFAST 
calculated  the  third.  That  approach  was  chosen  primarily  for  convenience  (FMASS  was  available 
1,0111  ,lie  data  and  AHC  was  known  from  the  type  of  fuel).  It  would  have  been  equally  valid, 
although  not  as  convenient,  to  use  either  ol  the  other  permutations. 

Regardless  of  the  manner  in  which  it  is  defined,  the  heat  release  rate  is  one  of  the  most  critical 
inputs  lor  CFAST.  Accordingly,  it  is  appropriate  that  we  consider  the  effects  of  changim*  this 
value.  Due  to  the  way  in  which  our  problem  was  originally  set  up.  the  easiest  way  to  effect  this 
change  was  to  adjust  the  mass  loss  rates  (FMASS)  and  allow  CFAST  to  calculate  new  values  for 


For  our  comparison,  the  Laundry  Room  model  was  run  with  one  half,  two  times  and  three  times 
the  nominal  heal  release  rales.  Figures  24A  and  24B  shows  the  effects  of  these  changes  on  the 
upper  and  lower  layer  air  temperatures,  respectively.  As  may  be  readily  seen. 'there  are 
significant  effects  on  the  qualitative  behavior  of  the  predictions  as  well  as  on  the  quantitative 
results.  \\  e  should  note  that  the  nominal  pyrolysis  rale  was  about  24  grams  per  second,  so  the 
results  shown  here  correspond  to  changes  in  pyrolysis  rate  of  only  a  feu  grams  per  second. 
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I'igure  24B.  Eilect  of  Heal  Release  Rate  on  Laundry  Room  Lower 
_ Layer  Temperatures _ 


5.1.2  The  OD  parameter 

products.  Clearly,  thiswHUffecT^ the  mediSte/soo^co8001  T  C.arbon  dioxide  in  ‘he  combustion 

will  also  affect  temperatures  (because  soot  is  an  important  factor  in ' SrCf  ^  °bvi°UsIy’  h 
concentrations  of  carbon  dioxide  and  rarhrm  r™  v,  ru™  tactor  ln  ra<hation  transport)  and 
unavailable  for  product ^  Wem  «  soo,  is 

the  Laundry  Room  model  using  OD  values  of  0  00^  06  and  0  m  Ii^  u  these  effects’  we  ran 
their  nominal  values,  as  given  in  Listing  T  °  '  ’  °'°6  and  °10'  AU  0ther  inPuts  were  set  to 

S“y  ROOm'  As  be  -peered. 

specify  fdSS  val^CFAST  wiU  predfcfnos''  ?f°  inpUt  3nd'  ‘h^^tf  the  use^fSk  to 
paragraphs,  this  is 

«r  ,*■ the  ^ 

which,  if  present,  cool  the  upper  layer  very  efficiently  The  w/8  >eiTfS?,V,ty  °f  S00t  Parlicles 
two  cases  in  which  soot  is  present is dZ  to  a the  temPera‘u^  for  the 
emitted  from  each  soot  particle  easily  escape  fro  i  l  ,  At  ‘°W  concent rations,  photons 

particle  density  becomes  high  enough  that  most  nh!  n  Y  bk’  a!,S^me  SOO‘  concentration,  the 
before  escaping,  increasingfhe 

concemrabons  "VS'™  ">c  soo, 

concentrations  in  the  lower  layer  lead  to  hthei  ruber  d  a,?  '°n'  M°re  ,mPortant|y>  high  soot 
to  the  fact  that  the  upper  layer  of  the  fire  a,r  ‘empera(ures.  This  is  due 

the  fire  plume  while  the  lower  layer  is  not  Asa  result'  d?'^  Y  h,eated.by  mass  injection  from 
•surroundings,  leading  to  a  net  loss  of  encro  v  hi  ,  ’  1 1C  'Plfr  layer  ,s  much  hotter  than  the 

d-  adjacent  upper  layer,  there  is  a  ,n  the  Iower  layCT’  which  is  cooler  than 

and  28.  Although  th^absollue  concentrations  of  th?  L,ailndry  Room  are  illustra(ed  in  Figure  27 
with  respect  to'oD  are  simi'a^K  T  Sp?C,esare  very  ^rem,  the  trends 

production  factors  are  seen  directly  the  concemVaiion  r'J  ?  e  ,ects  ol  a,tcnng  the  species 
with  OD.  As  mentioned  above,  this’ is  due  to  incre-isinc?™01 h  species  are  inversely  correlated 
soot,  leaving  less  carbon  available  for  any  other  species,' as  ODincreases°  "  °'C  f°rm  of 


:j  Note  .hat  tins  (approximate)  linearity  appl.es  onK  chan-es  in  the  OH  v  ,l„  .  ,  „ 

tZmi:' . . . "**“-* . . . 


model.  I  here  is  no 
regard  lo  pammeiers  oilier 
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Hgurc  26B.  Effect  of  OD  Parameter  on  Laundry  Room  Lower 
- -  Layer!  emperatures 
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For  the  lower  layer,  the  effects  of  changes  in  the  species  production  ratios  are  confounded  by 
other  factors.  In  particular,  lower  layer  concentrations  depend  on  mixing  effects  which  in 
CFAST,  occur  in  the  horizontal  vents  where  there  are  opposing  jets  of  hot  and  cold  air  entering 
and  exiting  the  compartment.  Since  vent  flows  are  affected  by  temperature,  pressure  and  layer 
heights,  the  lower  layer  concentrations  are  complex  functions. 

5.1.3  The  CO  parameter 

Figure  29  shows  that  carbon  monoxide  concentrations  are  approximately  linear  with  respect  to 
CO  in  the  same  way  that  soot  is  approximately  linear  with  OD.  Like  OD,  the  CO  input  has  a 
large  effect  only  of  the  species  which  it  directly  controls  and  relatively  small  indirect  effects  on 
concentrations  of  other  carbon-containing  species  (Figures  30  and  31).  However,  the  CO 
parameter  has  almost  no  effect  on  layer  temperatures,  as  shown  in  Figure  32.  The’practical 
implication  of  this  is  that,  unless  predicting  carbon  monoxide  concentration  is  an  important  goal 
of  the  simulation,  the  value  used  for  the  CO  input  is  rather  unimportant.  This  is  fortunate, 
because  CO  is  one  of  the  most  difficult  parameters  to  specify  a  priori. 

5.2  Geometry  Specification  Alternatives 

The  primary  reason  for  the  decision  to  use  a  compartment  rearranaement,  rather  than  a  virtual 
compartment,  approach  was  our  observation  that,  with  virtual  compartments,  the  Laundry 
Passageway  veilica!  vent  had  to  be  eliminated  from  the  model  in  order  for  the  simulation  to  run 
to  completion.  Using  the  rearrangement  method,  that  vent  could  be  included  (although  it  did 
have  to  be  reduced  in  size).  We  later  found  that  the  addition  of  the  Wardroom  lifted  this  vent  size 
restriction. 

In  this  section,  we  develop  the  "dummy"  compartment  concept  as  a  possible  method  for 
ciicum\enting  the  veitical  vent  pioblems.  We  also  illustrate  the  use  ol  virtual  compartments  and 
levisit  the  \eilical  vent  area  issue.  Although  not  all  oi  these  methods  were  used  in  our  case  study, 
we  believe  that  they  may  be  applicable  to  other  simulation  problems. 

5.2.1  Dummy  compartments 

Recall  that  the  virtual  compartment  method  involved  treating  a  sinsle  real  compartment  as  if  it 
were  two  tor  more)  "virtual"  compartments  (see  Figure  13).  The  "dummy"  compartment  concept 
differs  trom  this  in  that  we  proposed  to  add  to  the  model  a  new  compartment  that  had  no 
counterpart  in  reality. 

The  rationale  for  this  is  as  follows: 

1.  We  have  observed  that  stalling  problems  frequently  are  associated  with  the 
furthest  downstream  vertical  vent  (i.c.,  one  that  connects  to  the  exterior),  therefore 
it  is  possible  that  replacing  the  exterior  with  a  dummy  compartment  (so  that  the 
vent  in  question  no  longer  connects  to  the  exterior)  might  alleviate  the  problem. 

2.  We  know  that  the  effects  of  adding  a  new  compartment  propagate  upstream  (back 
toward  the  fire  compartment)  so  that  the  predictions  for  pre-exislimi 
compartments  can  be  altered  by  the  addition. 

We  hypothesized  that,  il  the  new  compartment  was  large  enough  and  had  a 
sill  I  iciently  low  flow  resistance  watli  respect  to  the  exterior,  the  new  compartment 
would  be  transparent  to  CTAS  1  the  model  predictions  would  be  the  same  as 
if  the  compartment  did  not  exist).  We  called  a  compartment  meet  inti  these  criteria 
a  "dummy"  compartment. 
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Figure  30B.  Effect  of  CO  Parameter  on  Laundry  Room  Lowei 
Layer  Soot  Concentrations 
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Passageway  and  addtria  neVcomp^^^  ^onntc^dto  version  °f  the  Laundry 

vent.  Since  we  had  no  guidance  S  6  LaUndr-V  Passageway  via  the  vertical 

compartment,  we  tried  10,  100  and  1000  mete?  (33  328  aphpropnaJe  Jor  a  dummy 

open  face  [areas  of  102  104  and  1  ne  m2  fflnnriv-  *  i.™  ,1  ^  ubes’  eacb  of  which  had  an 
comparison  with  our  prior  work  the  vertiLuenf  yJ°r  10*.and  10?  ft2)  respectively].  For 
companment)  was  set  to  an  area  of  0.25  nP  (2.69  fPU  “aHAra  mtde^86'™5’  ‘°  ^  dummy 

previous  results  in  HgureVI ^rseeThauh^pSdiied  i'165'  modelsfare  compared  with  the 

rxsxr  xkS  x,; 

predicts  air  temperatures  very  similar  to  those  obtained  without  the  dummy  compartmem.  nlet''°d 
oOhX^  c»e„  when  the  voiume 

the  fire.  This  suggests  that  the  concept ^mL  be  a  u  eh d -n  ^  ‘  S  ’  *"  lhe  of 

long  as  the  dummy  compartment  volume  is  not  excessiv^Cleariv  ?he'?eH  Circu‘f.lances  as 
some  inaccuracies,  especially  for  the  cornnarimernffm™'  a  ?,y’  techn>que  will  introduce 
those  errors  must  be  compared  with  tire  errors  ih ’ediately  upstream  of  the  dummy,  but 
vertical  vent  area.  P  'he  e,r0rs  lhal  would  be  caused  by  the  use  of  an  incorrect 

5.2.2  \  irtual  compartments 

-  >b=  combination  of  virtna, 


VT“',,U*  i equivalent  to  the  entire  8.51  m  (27  92  ft)  x  9  59  m  r«  so  n\  w  a  .  ,  ,  UUU1  10  me 

changes  permitted  the  Laundry  Passn«ew*iv  v^nirii  «.  .  ?  1  °  ^ ajdroom  bulkhead  .  These 

■«’  (8  40  fO)]  without  stallitm  die  ^mhuion  We  T  h  '"c;e;‘«d  to  the  correct  area  |0.78 
results  are  expected  to  be  different  than  those  fivu  S  °Mdi  note  l,iau  as  discussed  above,  these 

the  dummy  compartment  Unfortunately  due  in  t  h^W°U #d  have  beue.n  oblaincd  in  die  absence  of 
be  simulated.  unioitunately,  due  to  the  vent  area  problem,  the  latter  case  could  not 

«  might  expect  that  the  two 

true  of  the  lower  layers.  HoweS"? K .  £Ti?A<S  '?*  "Ould  he 

the  four  zones  independently  resulting  jn  ,i \Krnnt-  V  v  C  AS^  calcuIates  conditions  for  each  of 
compartments.  Tl.il  is  mZZK «  “o  C  boUnda?'.  boUVC“  «>™  v.rlnnl 
(figure  55A).  Note  that  the  two  values  br aekei  in'  .'  wpsi. nines,  which  diller  by  about  50°C 
Passageway  was  treated  as  Xml  «  R^lSr.'T  °^infd  "'"c"  "«  Lounrhy 
layers  (Figure  35B)  also  bracket  the  those  of  ih  *  •  i  °'  lL  'nU,a*  compartment  lower 

differences  among  (he  three  predictions  were  .^iigiblthuhLa^'1'  paSSageway  modd-  hul 
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Figure  33A.  Effect  of  Dummy  Compartment  on  Laundry 
Passagew  ay  Upper  Layer  Temperatures 

The  presence  of  a  dummy  compartment  changed  the  model 
predictions  for  the  Laundry  Passageway  upper  layer  but  there  was 
negligible  dependence  on  the  dummy  compartment  volume  over 
the  range  HP  to  1(P  nr  (35  x  10'  to  35  x  1(P  ft'). 
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Figuie  33B.  Ellect  o(  Dummy  Compartment  on  Laundry 
1  assageway  Lower  Layer  Temperatures 

Hie  piesence  oi  a  dummy  compartment  chanced  the  model 
■edm nous  (or  the  Laundry  Passageway  lower  layer  The  w as  no 

l()>  f  nC,|Ce,  n"  dummy  coniPartmeni  volume  below  !()<■  nr  (35  x 
there  were  ellects  at  the  largest  volume. 
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I  Jguic  34  B.  El  feet  of  Dummy  Compartment  on  Laundry  Room 
Lowe i  Layer  Temperatures 

I  he  presence  of  a  dummy  compartment  a  small  but  detectable 
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Eigure  35A.  Efiect  ol  Virtual  Comparlmenls  on  Laundry 
Passageway  Upper  Layer  Temperatures 

lempcraturcs  ol  the  two  virtual  compartments  that,  collectively, 
represented  the  Laundry  Passageway  were  different  from  each 
other  and  lrom  those  predicted  when  the  passagewav  was  treated 
as  a  simile  compartment. 
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Vertical  vent  size 


As  a  result  of  our  experiences,  it  became  obvious  that  some  compromises  would  often  be 
necessary  in  models  involving  vertical  vents.  We  have  shown  that  there  are  methods,  such  as 
dummy  compartments,  that  may  permit  correctly  sized  vents  to  be  included,  but  those  techniques 
introduce  their  own  errors.  Clearly,  the  user  will  have  to  make  a  trade-off  and,  to  do  that,  will 
need  some  idea  of  the  relative  magnitude  of  the  errors  of  the  different  techniques.  In  a  preceding 
section,  we  presented  information  regarding  the  inaccuracies  of  the  dummy  compartment 
approach:  in  this  section  we  consider  the  effects  of  altering  the  vertical  vent  area. 

The  impact  of  vent  area  was  investigated  by  systematically  varying  the  area  from  zero  to 
approximately  one  half  of  the  area  of  the  Laundry  Passageway  overhead.  Temperature 
predictions  for  selected  vent  areas  are  shown  in  Figures  36 A  and  36B.  For  the  upper  layer,  the  air 
temperatures  tended  to  rise  to  a  maximum  in  the  first  minutes,  then  abruptly  drop.  With  very 
small  and  very  large  vents,  a  quasi-equilibrium  state  was  reached  in  which  the  temperatures 
increased  slowly  for  the  remainder  of  the  simulation.  The  peak  and  quasi-equilibrium 
temperatures  decreased  as  the  vent  area  was  increased.  Venting  of  the  hottest  gases  via  the 
opening  in  the  overhead  would  be  expected  to  lower  the  air  temperature,  so  the  inverse 
dependence  of  upper  layer  temperature  on  vent  area  is  in  agreement  with  intuition. 

In  the  lower  layer,  the  situation  was  more  complicated.  For  the  smallest  vent,  the  temperatures 
simply  leveled  off  with  no  peak;  for  the  largest  vent,  there  was  a  pronounced  peak  and  a  large 
subsequent  drop.  The  data  are  inconclusive  for  the  intermediate  cases  due  to  the  stalls  which 
occurred  before  the  peak  temperature  was  reached,  but  suggest  that  the  temperatures  rose  faster 
(and  possibly  would  have  reached  higher  values)  for  larger  vents.  A  possible  explanation  for  this 
is  that  larger  overhead  vents  in  the  passageway  allow  a  greater  flow  of  hot  air  from  the  fire 
compartment. 

As  expected,  these  results  show  that  the  area  of  a  vent  in  the  overhead  has  a  significant  effect  on 
the  temperatures  in  the  compartment,  especially  in  the  upper  layer.  They  also  show  that  the 
magnitude  (and  even  the  sign)  of  this  effect  is  not  readily  predictable.  In  the  event  that  the  model 
does  not  run  with  the  actual  vent  areas,  it  may  be  necessary  to  substitute  a  different  area. 
However,  if  this  is  done,  it  is  important  that  the  user  carefully  investigate  the  effects  of  that 
change. 

6.0  RECOMMENDATIONS 

In  the  preceding  sections,  we  provided  examples  of  the  use  of  CFAST  in  modeling  a  specific 
problem  relevant  to  the  Navy.  In  this  section,  we  will  present  some  general  recommendations  for 
the  application  of  CFAST  in  other  problems. 

It  is  important  to  start  with  a  clear  idea  of  the  goals  of  the  modeling:  Are  you  trying  to  simulate  a 
specific  lire  or  are  you  interested  in  the  behavior  for  a  variety  of  possible  fires?  What  constitutes 
"good"  versus  "poor"  performance  tor  your  model?  How  important  are  temperature  predictions 
as  compared  with  species  concentration  predictions?  Is  the  interest  primarily  in  vertical  or 
horizontal  fire  spread? 

The  better,  and  more  narrowly,  the  problem  can  be  defined,  the  more  likely  it  is  that  useful 
results  will  be  obtained.  When  you  must  make  trade-offs  (and  you  wd]  be  forced  to  make  them  in 
all  but  the  most  trivial  problems),  refer  to  your  statement  of  the  goals  when  deciding  which 
details  may  be  neglected  and  which  must  be  kept. 
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further  discussed  here  because  examples  were  provided^ L  n^T "8  *  Spec,fic  fire  win  1101  be 
typical  case,  in  which  the  general  performance  under  a  variety  off-  °Ur  Study’  For  lIie  more 
use  of  a  Monte  Carlo  approach  is  suggested  With  th*t  ^  °^r^conditions  is  of  interest,  the 
various  ship  designs  could  be  compared8  u si n«  a  stadsbclunlw^  thfef^xPectfd  Performance  of 
the  results  of  that  analysis  could  be  presented  in  termJnfrn^K  °ufe  model  Predictions,  and 
given  a  specified  range  of  inputs.  For  examnle  mert^r  ^  Pr°babilities  of  various  outcomes, 
indicate  that,  for  fires  in  the  100  -  250  kW  range  there  is^.W^^h  ?.■“  beurthing  sPace  might 

entire  range  of  fire  sizes,  would  then  be  simnbmdS?  P  ^ A  senes  of  fires>  covering  the 

initially  be  described  in  general  terms-  Wha  are  the  flrTT  ofrlocations-  Possible  fires  should 
anticipated  fire  sizes?  In  what  locations  are  the  r  f  class,!, canons?  What  is  the  ranee  of 

possible  fire  types  or  locations,  it  is  probably  best  to  makeThem '°  °CCUr?  If  *here  are  multiple 
them  one  at  a  time.  F  y  St  t0  make  them  separate  problems  and  address 

ignored.  However,  i,  may  be  very  import™,  for  1  '1  Jl,sJraula,5  Wllh  CFAST  and  is  usually 
ventilation  effects.  If  so,  it  is  suooested  ilv.i  n-S  ?-P  b  and  you  ma>'  be  forced  to  include 
model  has  been  developed  and  successfully  run  T'hVnrim'1  SyStCni  be  added  after  the  rest  of  the 
that  development  of  the  mechanical  ventilation  mT>de|Pk  ,  reaSOn  °r  lhis  recommendation  is 
,l  “  P°in,leS-S  10  Peif°™  d-  work  until  tSSSrSS"  P,OCCSS  and 

time  to  start  buildin^die  actual  model"11 ^c^simp’le^c^*  b-een, accui?le|y  described,  it  is 
compartment;  therefore,  it  is  highlv  recommended  Iv  i  i  •  ‘f  lhaf  of  a  fire  in  a  single 
toward  accurately  defining  the  fire'and  the  «  modeling  effort  be  directed 

report,  the  lire  compartment  was  a  rectangula r  discussed  in  l<iis 

speed. cation  was  trivial.  In  cases  in  which  the  fire  ^  mn  f'.rc  compartment  geometry 

approach  is  to  start  with  a  simple  box  having  the  !  c  2m  ''!|l  *S  comP,ex'  a  reasonable 
concentrate  on  developing  the  most  realistic  f  e  !I!  r  lhe  actuaI  compartment,  and 

model  has  been  successfully  run,  adjul^!^  ^  th«  P^bminarv 

fue  compartment  geometry.  made  to  improve  the  accuracy  of  the 
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and  that  any  problems  encountered  by  addressed  immedi-i,HC  "n^  b<?  mn  a,,or  oach  addition 
"’C'cments  the  number  representin«  the  evi^  '  Y  -  Becausc  ,he  addition  of  a  new 
VVhM  and  CFCON  inputs  be  checked  to' ensure  S  m  U'°'’  I,n'X)rtant  l|uit  the  HVFNT. 
addition.  It  is  very  helpful  to  make  liberal  use  of  S,i"  COrrect  aflcr  each 
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Hie  area  specified  hy  a  V VHNt' bipmVand  Unu  "ilh  vcnt  sizcs  (especially  with 

wrdcr  10  pcrnm  ,hc  modd «« ™n-  ■"*  add->n  of  a 
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existing  restrictions  so,  after  a  new  compartment  has  been  added,  it  is  advisable  to  revisit  any 
previous  alterations  to  determine  whether  they  are  still  needed. 

When  the  correct  value  of  an  input  is  not  known  (for  example,  the  CO  parameter)  or  when  it  is 
necessary  to  substitute  an  incorrect  value  (as  was  done  with  the  WENT  area  in  our  test  case)  in 
order  for  the  model  to  run,  it  becomes  very  important  to  estimate  the  effects  of  this  inaccuracy. 
The  suggested  method  is  to  perform  a  sensitivity  analysis  by  running  the  simulation  with  a  range 
of  input  values.  If  possible,  these  inputs  should  bracket  the  entire  range  within  which  the  actual 
value  is  expected  to  lie. 

However,  in  some  cases,  such  an  analysis  is  not  feasible.  For  example,  when  we  initially 
attempted  to  apply  the  virtual  compartment  method  to  the  Laundry  Passageway  model,  we  were 
unable  to  find  a  non-zero  vertical  vent  area  which  did  not  cause  CFAST  to  stall.  Clearly,  it  was 
not  possible  to  conduct  a  sensitivity  analysis  under  these  conditions.  In  this  situation,  an 
appropriate  response  is  to  reconsider  the  overall  approach  to  the  problem.  In  our  example,’  we 
chose  to  use  the  compartment  rearrangement,  rather  than  the  virtual  compartment,  technique.’ 

Recall  that  CFAST  has  the  ability  to  reduce  the  heat  release  rate  in  the  event  that  it  calculates  that 
there  is  insufficient  oxygen  available.  For  this  reason,  it  is  useful  to  inspect  the  heat  release  rate 
output  to  verify  that  the  model  actually  used  the  values  that  were  specified.  This  need  not  be 
done  after  every  change  in  inputs,  but  should  be  done  (at  least)  after  the  model  is  complete.  Note 
that  deviations  from  the  specified  heat  release  rate  are  not  necessarily  incorrect  -  it  is  quite 
possible  to  specify  a  fire  size  that  cannot  be  supported  with  the  available  oxygen  supply. 
However,  deviations  should  be  investigated  to  rule  out  the  possibility  that  they  are  artifacts  of  the 
simulation. 


Time  =  1208.0  seconds 


Compartment 

Upper  Lower 

Inter. 

Pyrol 

Fire 

Pressure 

Ambient 

Temp.  Temp. 

Height 

Rate 

Size 

(Pa) 

Target 

(K)  (K) 

(m) 

(kg/s) 

(W) 

(W/mA2) 

I 

708.5  419.9 

0.50 

2.298E-02 

9.598 E+05 

-4.45 

1.156E+04 

2 

438.9  300.0 

1.2 

0.000E+00 

0.000E+00 

-0.849 

585. 

Outside 

O.OOOE+OO 

Table  9.  CFAST  Progress  Update 


A  C  PAST  progress  update,  similar  to  this  example,  is  displayed  on  the  console  device 
at  an  interval  (in  seconds)  which  specified  by  the  second  parameter  associated  with 
the  1 IMES  keyword.  II  the  F5  key  is  pressed  while  the  model  is  running,  the  current 
simulation  time  and  solver  step  size  is  appended  to  this  output. _ 

During  execution,  CFAST  displays  progress  information  (Table  9)  on  the  default  console  device 
(most  often,  the  computer  monitor)-'*  at  user-specified  intervals  (set  by  the  second  parameter  of 
the  I  IMES  keyword).  At  typical  CFAST  execution  speeds  (belter  than  10  time  faster  than  real 
time  is  common),  these  updates  scroll  very  rapidly.  However,  they  sometimes  slow  down 
drastically,  indicating  that  the  solver  has  made  a  bad  guess  regarding  the  step  size.  Normally,  this 
problem  is  quickly  corrected  and  CFAST  resumes  its  usual,  rapid  execution.  Our  experience  has 


-’'•  The  interval  between  progress  updates  (in  seconds)  is  set  by  the  second  parameter  of  the  TIMF.S  keyword.  Selims: 
this  parameter  to  zero  suppresses  updates:  the  default  interval  is  one  second 
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s  H  ?  indicate  * »« -nous 

are  indicative  of  a  stal^When  this  conditions0™™  operaPon-  Persistent  sub-microsecond  values 

edi,  the  input  fi,e  to  ^  Si™^ 

i“^o:S^tr,ali"g  ™ASS  vaiues  between 
t  e  fractional  opening  of  a  horizontal  vent  (CVENT  kevword'f  and^0  SOme  otber  inputs,  such  as 
However,  it  can  occasionally  be  the  cause  of  problems  ^n  narticnll™  f  ?i7  USefuI.ln  many  eases. 
(i.e.,  there  is  a  large  change  in  the  parameter  val.^  tw?  part‘cuIar,  if  the  transition  is  too  steep 

CfmST  may  encounter  numerical  problems  and  drasticall/slow^0  °Sely  Spaced  time  P°ints). 

p^dfS:.are  — *  *  * 

rrodU“  *  b*  -ing  RESTR  to 
work  a,  NRL  |3),  i,  was  desirable ,o  ^sSrnLu^ ZchT'*'10"'  examp^fin  previous 
fire)  to  a  mass  conflagration  period  (a  large  sprav  fireW  P-rebu.m  period  <a  sma»  Pan 

was  accomplished  by  simulating  the  pan  ffre  ndno  ,1  3t  3  specific  tIme  In  the  simulation.  This 
preburn.  The  model  was  then  ran  S w  S  t  S*  inpmS',f°r  lhe  d“™io„  of  tk 

force  the  simulation  to  resume  at  the  time  of  and  w  b?  S  “Sm?  RESTR  key'™'d  to 
the  first  run.  ’  dna  w,Ul  tbe  conditions  that  existed  at,  the  end  of 

obvious  question  w  ill  be:  How  good  are  the 
and  the  SHADWELL/688  test  results  In  between  the  model  predictions 

m  and  near  lhe  fire  compartment  (assuming  that  the  fi  e  w-i  rn^  ‘  If1  ag"eemem  ‘ends  to  be  best 
worse  as  the  distance  from  the  fire  compadinent  ncreases  bm  we  f3  defin<rd),and  Progressively 
which  the  predictions  for  one  layer  were  very  accnr  £ S  ,  f  f,'equently  have  seen  cases  in 
these  cases,  the  upper  layer  is  typically  better  than  lower  n  oT  b?f  1  ie  odier  layer  were  not.  In 
of  mixing  at  the  interface  between  layers.  ‘  ’  possibl>'  due  lo  the  absence,  in  CFAST. 

However,  it  is  important  to  understand  ilvu  cfact  n 

'oof  Realistically*  wc  can  only  expecl  CFAStC’  r'  a  models,  is  mi  ensinccrin.' 

cases:  d'is  is  snflicient  abso,“(e  valucC  In 

quesuons:  Is  ahernauve  A  likely  ,o  be  beuer  or  worse  than  B ».2S3ZS?""? 

7-°  SUMMARY 

WW  a  model  for  a  complex,  muhi- 
designers.  In  our  first  report  |6|.  we  showed 1  '0!'  ’  ’al  uould  be  available  to  ship 

hugely  on  known  or  estimated  fuel  and  combustio  n  proiirt!^  Onr''0'1  C°|Uld  b°  crca,cd  based 
the  geometrical  specifications  from  a  knowledge  of  t^  c  u't  ,  n  .?■  r"d  ,eport  |7<  developed 
we  have  provided  practical  guidance  to  persons  h  ,v LI  ,S  P  conf 'guralion.  In  this  report, 
fires  Due  to  the  inherent  limitations  of  the  CFAST  VoeaTnihrv^  i '°  3PP  y  CbAST  to  shipboard 
details  or  many  ship  compartments.  However  we  r  v  ,  y’,  'S  poss,ble  represent  the 

complex  lire  scenarios  using  reasonable  approximations.  '  ""  131  ”  poss,ble  to  describe  even 

from  hicran,rt!!ll,cs.h  I'K-raune  may  lx-  cstimak-d 

h;,st'd  or  ,lu- iiid mni:;;m“hy^:;“i,a,^\irl;t:lt™li^;; 
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parameters,  OD  and  CO,  cannot  be  accurately  estimated  from  a  priori  knowledge.  However,  we 
have  shown  that,  for  OD,  there  is  a  saturation  effect  and  that,  as  a  result,  many  fires  may  be 
treated  as  either  "clean"  (no  soot)  or  "dirty"  (sooty).  An  OD  value  0.06  appears  to  work  very  well 
for  the  latter  case,  which  is  the  one  of  primary  concern  in  Navy  problems.  We  also  demonstrated 
that  the  value  for  CO  is  often  unimportant  because  it  has  very  little  effect  on  temperature  or  on 
the  concentrations  of  carbon  dioxide,  soot  or  oxygen. 

The  techniques  used  for  modeling  of  complex  geometries  included  subdividing  and  rearranging 
compartments,  using  multiple  virtual  compartments  to  represent  a  single  real  compartment°and 
adding  dummy  compartments.  We  also  provided  examples  of  methods  for  estimating  mean 
values  of  thermophysical  properties  when  compartment  boundaries  are  composed  of  multiple 
patches  of  different  materials. 

We  found  that  some  situations  could  not  be  represented  at  all  (for  example,  it  is  not  possible  for 
one  compartment  to  simultaneously  conduct  heat  through  the  overhead  into  two  different 
compartments)  and  that  specific  combinations  of  input  parameters  could  cause  the  CFAST  to 
slow  to  the  point  of  being  effectively  unusable.  Because  there  is  no  known  way  to  predict  which 
input  combinations  will  cause  the  latter  problem  or  for  what  ranges  of  variables  this  problem  will 
occur,  we  suggested  a  methodology  for  resolving  stalling  problems.  This  approach  involved 
identification  of  the  offending  parameters  by  "commenting  out"  suspects,  adjustment  of  the 
parameter  values  until  the  model  runs  successfully  and  fine  tuning  of  the  parameters  to  be  as 
close  to  reality  as  possible  without  causing  CFAST  to  stall. 

It  has  been  found  to  be  very  advantageous  to  build  a  complex  model  by  starling  with  a  simple, 
working  case  and  adding  features,  one  at  a  time,  until  the  desired  scenario  is  reached.  When 
adding  compartments  to  a  model,  feedback  effects  can  significantly  alter  the  behavior  of  pre¬ 
existing  parts  of  the  model.  In  some  cases,  the  addition  of  a  new  compartment  may  permit  the 
elimination  of  approximations  that,  in  the  absence  of  the  new  compartment,  had  been  necessary. 
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Appendix  A 

CFAST  Keywords  Used  in  Modeling  of  the  Submarine  Ventilation  Doctrine  Configuration 


1 

VERSN 

Version  of  CFAST  for  which  the  input  file  is  intended.  A  title  may  also  be 
specified  for  identification  purposes. 

TIMES 

Simulation  time  and  frequency  of  on-screen  and  file  outputs.  ~ 1 

DUMPR 

Name  of  output  file  to  be  used. 

FTIME 

Event  timeline. 

RESTR 

Name  of  restart  fire,  if  any. 

Table  A-l.  Simulation  Control  Keywords 

, - - -  1 

IAMB 

Internal  ambient  temperature,  pressure  and  reference  elevation. 

EAMB 

External  ambient  temperature,  pressure  and  reference  elevation. 

WIND 

Wind  speed  versus  height  parameters. 
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Table  A-2.  Ambient  Environment  Keywords 
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DEPTH 

WIDTH 

HEIGH 

HI/F 

CEILI 

WALLS 

FLOOR 

THRMF 

HVENT 

WENT 

CVENT 

CFCON 


_  Depth  (x-dimension)  ot  the  compartments,  in  meters.  See  FPOS! 

-  Width  (y-dimension)  of  the  compartments'  in  meters.  See  FPOS.  ~ 

_  Height  (z-dimension)  of  the  compartments'  in  meters.  See  FPOS  " 

_  Floor  elevation,  relative  to  the  reference  elevation  in  - 1 - 

cdUngs.06  “  a"  entry  "'e  lhe™°PM“l  properties  database  describing  the 

Reference  lo  an  entry  in  (he  thermophysical  properties  dalabase  describing  the" 

"noorsenCe  10  a"  eMy  *"  ",e  lherm°p|'ysi«l  propenies  debase  describing  the  ' 

'  SS&S" te  U“d  'f  ^ified^IT- 

vem^ea  0110  S"lK  comPai1mcnts- lhe 

onhe  — 
— 


Table  A-3.  Model  Geometry  Keywords 
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LFBO 

Compartment  number  in  which  the  main  fire  is  located 

FPOS 

Coordinates  (right-handed,  Cartesian)  of  the  fire  location  within  the 
compartment  relative  to  the  lower,  left,  rear  corner. 

LFBT 

Fire  type.  Type  1  is  unconstrained  by  oxygen  availability;  type  2  is 
constrained. 

FQDOT 

Heat  release  rate  of  the  burning  fuel  at  the  times  specified  by  FTIME. 

FMASS 

Mass  loss  (pyrolysis)  rate  of  the  fuel  at  the  times  specified  by  FTIME  1 

CHEMI 

Miscellaneous  parameters  related  to  fuel  combustion  chemistry.  Includes 
molecular  weight  and  heat  of  combustion. 

HCR 

The  mass  ratio  of  hydrogen  to  carbon  in  the  fuel. 

02 

The  mass  ratio  of  available  oxygen  in  the  fuel  to  the  total  mass  of  fuel.  This 
applies  only  to  special  cases,  such  as  rocket  fuels,  in  which  the  fuel  consists  o 
a  mixture  of  oxidizing  and  reducing  agents. 

HCN 

Ratio  of  the  mass  ot  HCN  produced  by  pyrolysis  to  the  mass  of  fuel 
pyrolyzed. 

HCL 

Ratio  ol  the  mass  ot  HC1  produced  by  pyrolysis  to  the  mass  of  fuel  pyrolyzed. 

CT 

Ratio  ot  the  mass  ot  a  virtual  "total  toxics"  product  to  the  mass  of  fuel 
pyrolyzed.  This  product  is  taken  to  be  representative  of  the  combined  toxic 
elfects  ot  the  actual  pyrolysis  and  combustion  products. 

OD 

Mass  ratio  of  soot  to  carbon  dioxide  in  the  combustion  products.  This 
parameter  is  very  important  lor  correct  prediction  of  temperatures  (reference 

CO 

Mciss  ratio  oi  carbon  monoxide  to  carbon  dioxide  in  the  combustion  products. 

FHIGH 

Height  ol  the  base  ot  the  fire  (above  the  reference  position  established  by 

FPOS)  at  the  times  specified  by  FflME. 

FAREA 

Horizontal  area  ol  the  base  of  the  fire  at  the  times  specified  by  FTIME. 

CJET 

Switches  between  the  standard  and  ceiling  jet '  model  of  convective  heat 
transfer  to  the  ceiling. 
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Table  A-4.  Fire  Description  Keywords 

